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AB8T1ACT 

Th«  aec hen leal  In  (tract  Iona  between  ralnforcataanC  and  aatrlx  In  a 
coapoalte  undar  load  hava  baan  daacrlbad  by  strata  distribution 
aquations  in  tha  fibers  and  in  tha  natrix.  Although  apaclal  attantlon 
has  baan  given  to  strasaaa  resulting  fros  polyaar  last  ion  and  taaper- 
atura  shrinkage,  tha  solutions  obtalnad  can  ba  easily  adapted  for  usa 
with  ti'5  boundary  conditions  of  asternal  loads.  Several  cases  were 
considered: 

(1)  a  cylindrical  fllaaent  of  finite  length  aabedded 
in  a  raaln  cylinder; 

(2)  a  cylindrical  fllaaent  of  infinite  length  art  added 
in  a  raaln  cylinder;  and 

(3)  a  aatrlx  supporting  a  central  fiber  surrounded  by 
six  synaatrlcally  spaced  fibers,  each  in  turn 
surrounded  by  six  others. 


ill 


r 


F0B1W0BD 


This  final  report  was  prepared  by  Whittaker  Corporation,  Marne o 
Research  6  Development  Division,  San  Diego,  California,  under 
Contract  DA  44-177-AMC-20l(T)  entitled  "Mechanical  Kelationahip  of 
leinforcenenta  and  the  Binder  Matrix."  Tha  work  is  being  acconpliahed 
under  the  direction  of  B.  P.  McKinnon,  Contracting  Officer,  U.8.  Any 
Aviation  Materiel  laboratories  (USAAVLABS),  Port  Kuetle,  Virginia. 

This  report  covers  work  conducted  froa  10  June  1964  through  1  March 
1969. 

Work  on  thle  project  vaa  conducted  under  the  overall  direction  of 
Dr.  Juan  leaner.  Dr.  Gerhard  Wowak  contributed  as  a  consultant  to 
thla  program.  Report  has  been  rev  lowed  and  approved  by  B.  Levenets, 
Beeearch  Bagineering;  B.  Hldde,  Project  Office;  and  B.  Duft,  Manager, 
laglneerlng  Department. 


v 


ABSTRACT 


FOREWORD 

SYMOLS 


nmoDucnoM . 

TECHNICAL  DISCUSSION  . 

DtVRLOnCNT  OS  BASIC  EQUATIONS . 

BOUNDARY  CONDITIONS  OT  A  SINGLE  riBER 
Of  FINITE  LENGTH  IN  A  MATRIX . . 

SOLUTIONS  OT  THE  DIFFERENTIAL 
EQUATIONS  FOR  A  SINGLE  FIBER  IN 
THE  MATRIX . . 

APPLICATION  OF  TIB  SOLUTIONS  TO  THE 
BOUNDARY  CONDITIONS  EQUATIONS  FOR  TIB 
INTEGRATION  CONSTANTS . 

DETERMINATION  OF  THE  INTEGRATION 
CONSTANTS . 

THE  INFINITE  FIBER  IN  A  RESIN  CYLINDER  .  .  . 

THE  MULTIFIBER  PROBLEM . 

BOUNDARY  CONDITIONS  FOR  THE 
MULTIFIBER  PROBLEM . 

DIFFERENTIAL  EQUATIONS  FOR  DISTORTIONS 

SEPARATION  OP  THE  POTENTIALS . 

SEPARATION  OP  THE  VARIABLES . 

SOLUTION  FOR  (2  . 

SOLUTION  FOR  $3  .  . . 

COMPLEMENTARY  SOLUTIONS  FOR  % . 


CONTENTS  (Continued) 


!!C 

PARTICULAR  DfTECRAL  FOR  5  <C  ) .  108 

CITED  REFERENCES .  124 

UNCITID  RXPKRX1ICES .  125 

DISTRIBUTION .  126 

APPENDIX .  127 


viii 


SYMBOLS 


«o,  a  radlua  of  fibar  bafore  and  aftar  atraaa  davalopa  raapectlvely 

2aQ  radlua  of  ralnforcaaanta  bafora  atraav  davalopa 

a,  A,  conatanta 

b  dlanatar  of  raaln 

2bQ  radlua  of  raaln  cyllndar  bafora  atraaa  davalopa 

conatanta 

Blk  conatanta  daflnad  In  aqua t Iona  (182)  through  (217) 
conatanta 

D.  conatanta  daflnad  In  aqua t Iona  (218)  through  (223) 
and  aquatlona  (269)  through  (274) 

D..,  local  atraln  tanaor  Including  polyaarlaatlon  or 

J  tanparatura  contraction 

D..  voluaatrlc  changa,  or  dilatation,  or  traca  of 

tha  atraln  tanoor 


Dj,  atraln  tanaor,  axcludlng  polyaarlaatlon  or 

J  tanparatura  contract  Iona 

a1  unit  vac  tor  In  tha  dlractlon  of  partlcla  langth 

ap  unit  vactor  originating  at  point  f 

a11  unit  vactor  originating  at  point  Q 

K  Young 'a  uodulua  of  alaatlclty 

f^  daflnad  In  aquatlona  (224)  through  (268) 

df*  araa  alaaant 

g*  body  forca  danalty 

dG*  body  forca  danalty 

8ji  8*8*  tanaor 

H  Hankal  functlona 

k^,  kj  daflnad  In  aquatlona  (178)  and  (179) 


ix 


(Continued) 


2t0  length  of  specimen  before  etreee  develops 
itj»  U?"  length  efter  shrinkage 

sij,  sij  defined  In  equations  (180)  and  (181) 

n.,  n*  vector  nonsal  to  the  particle  area  In 
relnforceawnt  or  Matrix 

p^,  p  stress  vector  perpendicular  to  area  elaMent  df^ 


q  Index  referring  to  point  Q 
r  cylindrical  coordinate  In  the  radial  direction 
R  function  of  the  corrdlnate  r 
s  original  length  of  the  particle  observed 
?  length  of  the  particle  after  deformation 
dSj  component  of  dXj^  (u^)  along  *J1 


T 


Y  ur  u2,  u3 

V  V  V  U3 

X1 

xi  <V 


I  II 

a  #  a 

1 


8 


I 


or  8j» 

8. 


unit  vector  along  the  u.  coordinate  (tangential 

to  Uj)  J 

temperature 

general  coordinates  of  original  position  of  points 

general  coordinates  of  displaced  points 

location  vector  of  point  In  an  unstrained  SMdlun 

location  vector,  expressed  In  curvilinear  coordinates 
of  point  In  an  unstrained  medium 

principal  stress,  or  eigen  directions 

cylindrical  coordinates  In  axial  direction  (see  Figure  6) 

coefficients  of  thermal  expansion 

coefficient  of  shrinkage  due  to  reasons  other  than 
thermal  (e.g.,  polymerisation) 

coefficient  of  shrinkage 
eigenvalue  defined  In  equation  (430) 
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SYMBOLS  (Continued) 


•tJ-  «p 

X 

X 

Aijfct 

u 

v 

S1'  *i 
ij 

a  .  aAj 

all 

a12 

a13 

a21 

a22 

a23 

a31 

a32 

a33 
<P  t 

I 


defined  in  equation*  (225)  through  (233) 
Knonecker  delta 

atrain  tensor  due  to  expansion  or  contraction. 
Including  shrinkage  due  to  polymerisation 

a  Lam*  constant 

eigenvalue  of  the  wave  (defined  In  the  Appendix) 
tensor  of  elasticity 
.  / 

a  Lane  constant 
Poisson  ratio 

displacements  In  direction  physical  components 
of  the  dlsplaceaMnt  vector 

stress  tensor  component 
stress  In  radial  direction 
shear  in  tangential  direction 
shear  In  axial  direction 
shear  In  radial  direction 
atress  In  tangential  direction 
shear  In  axial  direction 
shear  In  radial  direction 
shear  In  tangential  direction 
stress  In  axial  direction 
cylindrical  coordinate 
potential  function 


SYMBOLS  (Continued) 


Superscripts : 

I  refers  to  reinforcement 

IX  refers  to  matrix 


Subscripts : 

l»J»k*4,r,  or  s  free  index 


s  set  of  orthogonal  unit  vectors 

g  metric  vectors  tangents  to  the  curvilinear 

^  coordinate  elements 

dUj  curvilinear  coordinate  elements  Figure  (l.b) 

metric  tensor  of  the  strained  body 

ds  differential  arc  length 

V  displacement  vector 
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This  program  was  initiated  to  provide  further  understanding  of  the 
mechanical  Interaction  between  reinforcements  and  the  matrix  in  a 
fibrous,  composite  material.  The  mathematical  fundamentals  Initially 
established  were  a  set  of  solutions  of  partial  differential  equation 
describing  the  distortions  In  both  the  reinforcement  and  matrix.  The 
solutions  were  kept  In  general  form  so  that  they  could  be  used  In 
establishing  those  internal  stress  distributions  resulting  from 
polymerisation  and  temperature  shrinkages  as  well  as  those  resulting 
from  loads  Imposed  at  the  outer  surface. 

Two  jases  are  considered:  (1)  a  cylindrical  filament  of  finite  length 
centrally  embedded  In  a  matrix  cylinder,  and  (2)  a  matrix  which  supports 
a  central  fiber  that  Is  surrounded  by  six  symmetrically  spaced  flbere. 

Undulating  stresses  apparently  result  In  both  cases.  Complete  axial 
symmetry  was  assumed  In  the  monofilament  case,  while  for  the  seven- 
fiber  model,  the  stress  pattern  Is  repeating  six  times  around  the 
central  fiber. 

The  belief  is  that  the  first  case  has  been  solved  rigorously  In  this 
work  without  any  simplifying  approximations  or  assumptions.  For  the 
latter  case,  certain  potential  functions  representing  the  component  of 
the  curl  of  the  displacement  vector  and  one  function  of  volumetric 
compressibility  have  been  set  to  sero;  to  date,  however,  there  Is  no 
rigorous  theoretical  support  for  this  assumption. 


INTRODUCTIOH 


The  strength  of  composites  appears  to  ba  wall  balow  that  which  might 
ba  achieved,  when  the  available  strength  of  the  reinforcing  filaments 
la  coneldered.*  Any  attempt  to  develop  reinforced  composites  with 
etrangth-to-welght-ratloe  higher  than  those  presently  attainable  must 
be  baaed  on  store  complete  Information  concerning  the  mechanical 
Interaction  between  the  reinforcement  and  matrix  In  a  composite.  The 
subsequent  text  susssarlaes  the  mathematical  development  accomplished 
under  a  program  whose  goal  was  to  ss>re  clearly  define  this  Interaction. 

The  Initial  consideration  la  that  of  distortions  In  two  hoswgeneous 
materials,  one  esfcedded  In  the  other.  A  consideration  of  equilibrium, 
a  stress-strain  relation,  and  a  straln-dlsplacesmnt  relation  will  lead 
to  differential  equations  for  the  stresses,  strains,  and  dlsplacesmnts 
on  and  within  a  solid  body.  Aa  a  specific  caae,  a  cylindrical  system 
la  considered  In  soon  detail.  This  basic  development  Is  then  applied 
to  two  cases. 

Case  1  -  A  cylindrical  fllaamnt  centrally  embedded  In  a  matrix  cylinder 
la  studied.  Axial  symmetry  Is  assusMd.  The  sMtheswtlcally  rigorous 
solution  of  the  three  differential  equations  of  displacement  Is 
derived.  The  proper  boundary  conditions  applied  to  the  solutions 
results  In  the  equation  for  strain  and  stress  distributions  In  the 
composite.  Special  attention  Is  given  to  the  boundary  conditions  on 
the  Interface  between  reinforcement  and  matrix,  and  to  the  case  where 
external  forces  are  applied.  In  this  manner  the  residual  stress 
distribution  In  the  composite  Is  obtained.  A  specific  numerical 
exaaple  of  an  Infinite  fiber  surrounded  by  resin  Is  given. 
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Case  2  -  A  matrix  supporting  •  central  fiber  that  la  surrounded  by  six 
aymamtrlcally  spaced  fibers  Is  treated.  The  mathematical  steps  are 
slaillar  to  thore  followed  In  the  first  case,  but  ere  much  more 
c oeq>llcated  as  there  Is  only  hexagonal  syammtxy  around  the  fibers.  This 
and  other  geometric  limitations  are  assisted  as  boundary  conditions.  The 
boundary  conditions  on  the  Interface  between  the  relnforceamnt  and  the 
compoelte  are  the  seam  as  those  In  the  single* fiber  case.  The  basic 
equilibrium  equations  (expressed  In  terms  of  displacements)  are  con* 
sldered  simultaneously,  and  general  solutions  for  the  displacements  or 
distortions  In  the  composite  ere  obtained. 
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The  coaponent  aster la Is  of  the  composite  are  anuMd  to  bo  undor  intornol 
and  external  atrosaos.  Included  art  polymerisation  shrinkages  and 
temperature  expansion  or  contraction  phenomena.  The  total  strain 
resulting  from  such  stresses  Is  represented  by  the  strain  tensor. 

In  cylindrical  coordinates,  the  eleaents  of  the  strain  tensor  are 
obtained. 


*  Bquatlon  (A71)  (A72)  In  the  Appendix  to  this  report. 
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Equations  (1)  through  (5)  introduced  into  ths  stress  tensor* ** 


result  In  stress  distribution  In  radial,  tangential,  and  axial  direction, 

As  symbolised  In  the  following  o..^  ,  ,  o^  direct  stresses  In 

radial,  tangential,  and  axial  direction,  and  »  ^3  »  th* 

shear  strasses  In  tangential,  axial,  and  radial  direction. 

(Mote:  All  of  the  stresses  are  written  with  subscripts,  which  Is  custom¬ 
ary  for  the  theory  of  elasticity.  Superscripts  Z  and  ZZ  Indicate  the 
stress  In  the  relaforcesmnts  and  resin  respectively.) 

to  obtain  the  displacement  f  J  .  ,  |*  and  «“  .  |“  , 

the  differential  equations  of  displacements  must  be  solved.  Restated, 
the  partial  differential  equations  are 

Zn  the  radial  direction: 


*  Equation  (63)  In  the  Appendix 

**  Equations  (73)  In  the  Appendix 
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In  the  tang«nt 1«1  direction: 


a2*j  !  !  a2?! 

,2  +  2(1  -  2v)  r  at  a<*>  + 

o* 


2(1  -  2v)  r  dCP  ds 


♦ 


J- 

2r 


l  3  -  4v  25i+  !ii.  l 

r  1  •  2v  W  dr  r 


(aa$+w}  +  {«!  •ln'f +  *2 co,!()  H*  *  0  <9> 

In  the  axial  direction: 

1  i!*i  j_  i!h  i-v  a2**  i  a2<i 

2  at2  +  2t2  a*2  +  1  •  »  a.2  *  2<1  *  w)  ar  a*  + 

i  .  *2  ,i  i  a«i.  i_  aJa  . 

2(1  -  2v)  r  dtp  da  2(1  -  2v)  r  ds  2r  dr 


m  t>K  H)  h*  • »  »•> 

It  can  be  aeen  that  the  three  differential  equation*  are  coupled  partial 
differential  equation*  which,  by  Introduction  of  adequate  aatheeatlcal 
potential,  can  be  aeparated.  Before  this  1*  considered,  however,  an 
extensive  study  of  the  boundary  conditions  Is  presented  to  provide 
Insight  to  the  problem  and  facilitate  its  solutions.  The  first  Is  a 
finite  fiber  In  a  astrlx  cylinder,  and  the  second  condition  consists  of 
a  matrix  and  a  cent  al  fiber  surrounded  by  six  symmetrically  placed  fibers. 
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BOUNDARY  CONDITIONS  OF  A  SINGLE  FIBER  OF  FINITE  LENGTH  IN  A  HAHU1 


Consider  a  fiber  of  diameter  2a  and  length  2<t  embedded  In  a  matrix 

o  o 

cylinder  of  outer  diameter  2b  and  length  24  .  The  Index  tero  refers 

o  o 

to  dimensions  before  any  stress  has  developed. 

♦. 


Figure  1 


The  systeai  of  the  partial  differential  equations  (8)  through  (10) 
describes  the  distortion  in  both  materials.  It  is 


«J(r,  «) 


(11) 


the  distortion  vector  In  the  fiber  and 

{“(r.  .)  (12) 

the  distortion  vector  in  the  matrix.  The  material  constants  E  ,  v  , 

8  ,  a  *  which  have  superscripts  I  or  11  refer  to  the  reinforcement  or 
matrix  respectively. 
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The  fundamntal  equations  of  the  theory  of  elasticity  are  valid  and 
the  results  obtained  in  the  Appendix  can  be  utilised.  The  boundary 
conditions  ere  laplamnted  by  the  fact  that  both  enda  are  free 

froai  stresses. 

If  there  are  no  external  forces,  then 

°ij(bu  ’  *)"j 

is  considered  a  set  of  boundary  conditions. 

•!!(“.■  ■  1 

is  the  etreas  tensor  (equation  7)  on  the  aurface  of  the  aatrlx  in  the 
points 

(r  *  bn  ■  b.(l  - !  *) 

and  nj  la  the  unit  vector  norm l  to  the  aurface  elemnt  in  the 
point  bj  ,  s 

The  length  bjj  ■  bQ|  1  -  |^J,  the  outer  radius  shortened  by  the 
shrinkage  of  the  aatrlx.  But  alnce  n^  la  a lvays  in  1  (or  r)  direc¬ 
tion,  It  Is,  in  the  aesuaed  cylindrical  coordinate  syetea, 

°J 

so  that  equation  (13)  become 

(bn  • fllJ  6u 


0  (13) 

In  equation  (13), 
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or,  more  specif ically, 


°ii  |bn  •  *)  *  0  <“> 

°\\  (bn  •  •)  *  0  <l7> 

r31  |bII  •  *)  ’  0  <l4> 


By  using  the  equation  for  the  stress  tensor,  these  boundary  conditions 
result  in 
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_b_ 

dr 


0 


(23) 


Equations  (22)  through  (24)  can  alto  bt  used  for  the  Infinite  fiber. 
The  boundary  conditions  for  the  fiber  at  the  end  of  the  specimen  are 


°ls  (r  •  "j  ‘  0 

°li  lr  •-<!)  “j  '  0 

and,  for  the  aatrlx,  are 


for  0  <  r  i  a  (25) 


*“('  •  •  0 
•t )  lr  •  <ti  “J  ‘  0 


for  a  *  r  *  b 


(26) 


In  equations  (25)  and  (26),  consideration  oust  be  given  to  the  fact  that 
the  length  K  In  both  cases  Is  subjected  to  shrinkage,  so  that 

«l  ■  *0  I1  • 

and 
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Hare,  th«  noml  vector  la  In  a  direction  (exiel) 


"j  *  *3)  ‘  *3) 

Equetiona  (25)  end  (26)  reault  in  r he  conditiona  for  the  reinforceattnt: 

°33  { r  '  ®tj  *  *  °33  (r  •  *  ‘  0  (2,) 

'It  lr  •  <i)  ‘  *  'It  (r  •  *  <i)  ■  0  «•> 

a31  (r  *  ^l)  *  *  03i  (r  »  *  *4)  “  0  (29) 

end  reaulta  in  the  following  conditiona  for  the  aetrix: 

°33  (r  1  ^ll)  -  -  c\\  (r  .  -  tjj)  •  0  (30) 

c32  ('  •  <ti)  -  *  'll  («  •  -  Si)  •  0  <3l> 

'll  I*  •  hi\  *  *  °S  lr  •  *  <nl  ■  0  <32> 

Applying  equation  (7)  again  yielda 

D33  (r  •  tXll  ♦  T^T  Dtt  lr  •  *  0  <33> 

®32  (r  ’  4<l)  *  0  «*> 

D3I  (r  *  14l)  ‘  0  <«> 
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DJ3  ( r  '  **11}  +  l  .VJVII  Btt  |r  •  **11) 

D33  I*  •  ‘<11) 


ud  with  equations  (1)  through  (6), 


S1  I*  •  “l|  *  TT1-  I'  •  “J 


£  I'  ■  “-I} 


(37) 


Dl\  I*  •  ‘Si)  •  0  <M> 


^2  I  1 

dt  r  •  ±4tiJ 


dr  (r  •  *  ds  (r  *  "  0 


*  r^(}  ^  I- •  «,,!* 

?r  |r ,  *<„)) . 
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■  0 


(43) 


ds 


(r  •  ±<ll) 


dr 


(r  *  ±<ll) 


-  0 


(44) 


Equations  (33)  through  (44)  are  twelve  boundary  condition*  relating  only 
to  the  finite-length  fiber.  These  conditions  mist  be  sat.sfied  by  the 
solutions  of  the  differential  equations.  To  define  the  problem  for  the 
infinite  fiber,  the  length  l  mist  be  assumed  to  approach  infinite, 
which  is  a  separate  astheaatlcal  problem. 


For  reasons  of  symmetry,  the  following  conditions  are  valid: 

I  I 

5“  (r  .  .)  •  (r  .  -  «) 

I  I 

Sj1  (r  ,  s)  -  -  (r  .  -  s)  (45) 

1  I 

-II  ,  ,  -II  ,  . 

*3  (r  ,  s)  •  ’  53  (r  ,  -  «) 

Since  in  the  plane  s  ■  Q  ,  equation  (45)  must  be  satisfied, 
can  be  developed  and  represented  by  Taylor's  series. 


*) 


0) 


JL 

a." 


0) 
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or 


0)  -  0 


(46) 


Proa  equation  (46),  it  follows  that  for  1  -  1  ,  all  uneven  partial 

derivatives  with  respect  to  s  at  the  location  (r  ,  0)  must  vanish 
Identically.  All  even  partial  derivatives  with  respect  to  z  vanish 
for  1  ■  2  and  3  ,  so  that 


and 


(r 


0)  -  0 


I, 

2,3 


(r  ,  0)  -  0 


(47) 


The  solutions  of  equations  (8)  through  (10)  saist  satisfy  equation  (47). 

A  third  set  of  conditions  can  be  derived  from  t.  .  face  that  on  the 
fiber  matrix  Interface,  two  neighboring  particles  of  the  two  materials 
aust  reset  in  together  during  the  distortion  displacements.  First, 
assume  that  both  the  fiber  and  the  eatrlx  are  made  of  an  homogeneous 
am  ter la 1.  Secondly,  assume  that  the  shrinkage  of  a  cubic  particle  of 
the  length  1  Is  1-0*  and  1  -  0**  respectively.  Third,  assusm 
that  0**  >  0*  .  Then  the  fiber  will  also  be  compressed  in  all  direc¬ 
tions.  As  a  result,  stresses  will  develop  In  both  siaterlals.  Before 

the  stresses  develop,  the  common  dimensions  are  a  ,  b  ,  and  l 

o  o  o 

If  the  fiber  were  not  present,  the  Inner  radius  of  the  resin  would 
shrink  freely  until  It  would  assume  the  dimension 


(48) 
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Also,  the  outer  radius  of  the  resin  reduces 


b 

o 


I 


1  -  B 


and  the  length  be  com  s 


(49) 


(50) 


A  particle  in  the  matrix  with  the  original  coordinates 


Similar  results  are  obtained  for  the  coordinates  of  the  reinforceMnt : 

a1  -  ao  j 1  -  B1 )  (51) 

l1  -  (l  -  81)  (52) 

A  glass  cylinder  in  a  hollow  resin  cylinder  is  depicted  in  Figure  2  in 

order  to  illustrate  the  distortions  caused  by  the  difference  in  shrink- 

age  and  thermal  expansion  of  the  two  materials.  Each  material  by  itself 

would  shrink  to  the  cylinder  sice  indicated  by  the  thin  lines,  but  is 

restricted  in  this  case  ly  the  presence  of  the  other  material.  Since 

material  particles  must  stay  in  contact  in  the  Interface,  the  point 

P„  (in  the  state  without  stress)  settle  st  the  common  location  P 
o 

Instead  of  moving  to  the  locations  (for  the  reinforceMnt)  and 

Pjj  (for  the  matrix). 
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(a)  Diagram  Depleting  Compo 
nanta  before  and  after 
Independent  Shrinkage 


(b)  Composite  Shape 
Resulting  from 
Combined  Shrin¬ 
kage 


Matrix  before 
Shrinkage 

Compel  1 t» 
Interface  before 
Shr Inkagc 

Se Infer cement 
before  Contrac¬ 
tion 


Free  Reinforcement 
after  Contraction 


Composite  Interface  ifter 
Shrinkage 

deformed  Composite 


(c)  View  A  (Enlarged)  Showing  Combined  Shrinkage 
in  Relation  to  Independent  Shrinkage  of  Both 
Materials 


Figure  2.  Reinforcement  Cylinder  Embedded  in 
Hollow  Matrix  Cylinder 
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f 


The  total  distortion 


can  ba  described  in  tv^  ways: 


(53) 


and 


P  P  -  P  PT.  ♦  PTT  P 
o  o  II  II 


(54) 


Pg  P  is  the  distortion  of  the  reinforceamnt  cylinder  as  in  the  points 

Fi (*•  I1 '  #II  • 


There fora, 


■  •I1-®1! 


Similarly,  distortion  for  the  matrix  is 


'll'  •  (  *o  I1  -  ®IX)  •*l1-'11) 


(35) 


(56) 


The  vactors  P  P.  and  P  P__  are  displacamsnts  of  P  a  ,  a  in 
O  1  O  11  o  o 

the  interface  of  Z  and  II,  so  that 


F  P, 

o  I 


'll  •/  ’  *3i  ,pI 


(57) 


and 


ft  •  V11 
o  rII  vi 


"  *li  V"  ’  *3i  ,plX 


(56) 
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Using  equations  (S3)  and  (54),  and  equations  (5S)  through  (58),  the 
following  relation  exists: 


or 


<  [•„(*  *  bI)  •  ‘I*  - sI)]*  5“[*oll  •  8“)  >  *('  - e“) 

*  *11  ‘el*1  ‘  8“)  +  *31  'I8*  *  S'1)  <”> 

The  subsequent  boundary  condition  equations  follow  from  equation  (59): 

<$.|>  •  *1  •  ■!*  ■  >'j  •  *“[-.(>  •  *“l  ■  -I'  •  -'l] 

-  •  I*1  -  011)  (60) 


4.i‘  •  *i  •  +  •  ‘‘i]  •  •  ’"i  •  -i'  ■  •“)]  •  • 

>;[■.('  •  .‘I  •  •  *‘l]  •  >!'[•.['  •  •“] .  •(.  -  .“lj 


(•■  •  *i 


Using  the  transformations 


r1  -  r  (l  -  31) 


(61) 


(62) 


(63) 
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t1  -  r  (l  -  B1)  (64) 

In  material  I,  and  th«  transformations 


r11  -  r(l  -  B11) 

(65) 

-  »(l-  811) 

(66) 

In  material  II,  equations  (60)  and  ('  then  becoaw 

l\[.1  .  .*)  -  -  .» 

(67) 

<‘|.1  .  .*)  -  .  .“)  -  .»  - 

(68) 

By  taking  into  account  a  strip  PQ  and  its  elongation  or 

during  the  process,  and  considering  Figures  1  and  2(c). 

contraction 

Qx  ♦  JFJ  ♦  5^5 

(69) 

But  it  is 

"T  -  *  -  »x|  •  •(>  - 

and 

v  *  •  bI)  •  (* +  HI1  *  #Ij  ■ 
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Writing  If  ■  V*  ,  then  equation  (70)  becomes 

vj  -  a«(i  -  ^  ^  (.o(i  -  e1)  .  ,(i  -  b1))}  (71) 

and  for  the  matrix  becomes 

<?•<■(■  -  *  'S  (•.!'  •  -“I  •  •  ■"§  <’» 

Since 

vj  -  vj1  (73) 

it  follows  that 


Equation  (74)  could  also  have  been  obtained  by  differentiation  from 
equations  (60)  and  (61);  this  method  can  be  used  to  verify  the  procedure. 
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f 


I 


On  Che  interface  of  fiber  and  retin,  Chere  are  Che  normal  acreasea  in 
both  mate rial a  equal  and  oppoaiCe 


and 


•  ■(>  ■  !,l)*u  >”> 

•  •  *“i .  -i>  •  •■■))*„  <»> 


(77) 


Consequently,  the  stress  boundary  condition  at  the  interface  is 

-:,(•>  •  ••)  ■  -  «'i)-  -  *“i  ■  *(■  •  *■'(  •  •  <»> 

expressed  in  coordinates,  equation  (78)  is 

‘  S‘l  •  *(*  *  S1)]-  -  BU)  ,  l|l  •  I11)]  -  0  (SO) 

»!,[-.(■  -  »‘l  ■  -I*  •  •*)]-  •!![•.!■  •  >“) .  • »“)]  ■  •  i«» 

The  six  boundary  conditions  on  the  common  surface  of  the  fiber  and  resin 
are  in  equations  (60)  through  (62)  and  equations  (79)  through  (81). 
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Solution  of  Hi  Differential  Equations  for  j  jj^j  Fiber  la  the  Matrix 

After  having  attabliahad  all  boundary  conditions,  tha  differential 
equations  (9)  through  (11)  are  treated.  The  differential  equations 
for  one  single  fiber  are  brought  into  the  forne 


Lai  +  _ai  + 1  .a.)  ( i!*i  +  i!*a  + 1  !S 

la.2  ar2  r  *rIW  »»2  r8r 


(•A) 


Vlret,  equation  (82)  is  solved  by  setting 


2>2«i  a2«, 


ds2  dr2 


*  Mr-‘) 


(85) 


wh.r.  0  wilt  b«  .  .olutlon  of 


i!*i  ♦  2jh  *  i  2*1 .  X  # 

a.2  ar2  * r2  *l 


(86) 


According  to  iernolll,  this  equation  can  be  solved  by  setting 


-  A^(r)  ^(s) 
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*r 


and  the  following  ordinary  differential  equations  obtained: 


0 


with  the  solutions 


sad 


Hence,  equation  (63)  la 


E^r)  ^(s) 


(87) 


(88) 


(89) 


(90) 


(91) 


Using  a  procedure  slallar  to  that  used  for  equation  (86),  the  'allowing 

Is  set 


-  E(r)  Z(s) 


(«2) 
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And  obtained  (Ref.  2) 


10^1 


ds 


Two  sepereble  ceeee  exlee  w^en 

Z1(«)  -  a  Z(«) 


(93) 


(94) 


end  when 


Rx(r)  -  |  R(r) 


(95) 


Both  canes  ere  possible  solutions  end  the  sub  of  both  is  another 
solution.  For  equation  (94) t  it  follows  fron  equation  (93) 


1  ^ 


lfdit 

»V  rdr 


-  “V  R  - 
2 
r 


aRjOr)!  • 


(94) 


fron  which  It  follows  that 


♦  kj  Z  -  0  (97) 

ds* 


and 


d?L  +  l  ik 

dr2 


+ 


m  -i-  R 

*1  2  * 
r 


a  R|(r) 


(98) 
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Because  of  (67),  (94),  and  (97),  tha  following  axlats: 


k 


1 


-  \ 


(99) 


and  equation  (98)  it 


.  2  r  dr 
dr 


-  X2  -  4  R 

i  r 


a  ^(r) 


(100) 


g^(r)  In  equation  (100)  la  a  linear  combination  of  the  aolutiona, 
defined  by  equation  (90)  aa  Bessel  functions.  The  total  solution  of 
differential  equation  (100)  is: 


R(r)  -  J^(iX.r)  ♦  H^UXr)  ♦  C^r  Jo(Ur)  +  D^r  Ho(iXr) 


(101) 


For  the  first  case,  or  equation  (94),  the  general  solutions  are 
*1  *  [A1  Ji<1X,)  +  B1  H{l)<1Xr)  ♦  V  J0(Ur)  ♦ 

V  Ho(lXr)j  [.1X*  ♦  .*1X*j 


(102) 


For  tha  second  case,  it  follows  from  equations  (95)  and  (93)  that 


l  d2z  zi<*> 

*H-» -4- 

ds 


-  K 


(103) 


ii  +  ia*  |k2 .  i.| ,  .  o 

dr2  P  r2 


(104) 
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Because  of  equation  (95),  the  equations  (104)  and  (SB)  are  identical,  thus 


k2  -IX 


(105) 


and  equation  (103)  becostts 


di 


a  «U*  ♦  8  .'U* 


(106) 


Two  Independent  solutions  of  equation  (106)  are  already  found  in 
equation  (89).  The  two  other  solutions  can  be  found  by  setting 


Z(«) 


Y  u(s)  eiX‘  ♦  6  v(z)  e’iX* 


(107) 


Equation  (106)  with  equation  (107)  becones 


y(u11  ♦  2  uli\  -  X2uJe^X*  ♦  fijv11  -  2  v^X  -  X2v]e 


-IXs 


a  etX‘  ♦  8  e‘lX*  (108) 


From  equation  (108),  the  following  Is  obtained: 


11  ,,  i.. 

u  +  2  u  IX 


11  , 

v  -  2  v  IX 


.2  ,  a 

Xu  u  - 


XJv  ■  £ 


(10S) 

(110) 


Since  a  end  0  ere  cos^lex,  It  follows  that  for  the  laaglnary  part 


ul  -  *  i«ag  ( ?  I  ■  C,  A 


2X 


2  Y 
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1 


Consequently, 


^2 

v(«)  -  -f  t 

The  reel  pert  gives  solution*  elreedy  known,  end  equetion  (107) 

becoM* 


(111) 


(112) 


2(e)  •  Cjte*1**  ♦  d2*e"1X*  (113) 

The  complete  solution  of  the  ordinery  differentiel  equetion  (106)  ie 

Z(e)  -  e2  e1X*  ♦  b2  e~1X*  ♦  C2*  eiX*  ♦  d2*  e“1X*  (114) 

end  therefore  the  totel  solution  for  the  second  cess 

•  Ja2  J^(iXr)  ♦  B2  H<l)(Ur)J  £elX*  ♦  e2  e’iX‘  ♦ 

b2*  elXr  ♦  C2*  elXrJ  (115) 

The  totel  solution  including  both  ceses  is 

?!  -  ?!  +  ?!  <1U) 
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or 


MH 


|l2  J^(lXr)  *»-l2  H<l)(lXr)j£ 


#u.  + 12  ,-U.l  ,  + 


^  J^(lXr)  +IL  H<l)(lXr)  ♦  J0(lXr)  ♦ 

H<l)(lXr)j  •  j«lX*  ♦  il  e‘iX*j  (117) 

Using  th**  mm  procedure  for  the  distortion  vector  conponent  In  exlsl 
direction,  differential  equation  (84),  a  solution  analog  to  equation  (117) 
Is  obtalued: 

f3|r,s)  -  ^i2  Jo(l?r)  +  I2  H^l)(lfr)J^* -►  I2  e‘l^J  s  *»■ 

^  Je(llr)  ♦  Hj^dSr)  ♦  8jT  JjdXr)  ♦ 

Bjr  H{l)(lltr)J  Je1^*  ♦  a1  e'1^  (118) 

The  general  solution  of  equation  (83)  Is  similarly  obtained;  this  Is 
a  partial  differential  equation  which  results  In  a  Bessel  equation  of 
first  order  and  a  differential  equation  of  exponential  functions 


4*±a  +  |2  xU  .  o 


72T7d7■,, \»  2 

dr  \  r 
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and 


*V*2z  •  o 

d*Z 

The  solution  of  Aquation  (83)  is  the  sum  of  the  particular  solutions; 

i  •  A  •  ^ 

■  (v +  t1]!1  +  *i*) +  (*2  vtJr> + 

B2  H^l)  (iJr)J  |a1^*  ♦  I2  .-*)  (119) 

The  solutions  for  the  distortion  vector  components  in  the  reinforcement 
(superscript  I)  and  the  matrix  (superscript  II)  are  aa  follows: 

In  Radial  Direction  in  the  Reinforcement: 


$J|r,*J  ■  J1|iXIrJ  |aJ  cos  XXs  ♦  aJ  sin  XX«J 


s  ♦ 


rl(iXXr)  ♦  aXr  Jo|iXXr)j  JaX  cos  XXs  ♦  aJ  sin  X1^ 


(120) 


In  Radial  Direction  in  the  Matrix: 


*  ^(u11.)  +  .«H<l>(«IIr))(A“ 


cos  XXX«  ♦ 
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♦  a“  sin  XIXe|  s  ♦  | JtXIXrJ  ♦  eXI  ll[l*|iXX1r)  + 


•?*  JollxlIr)  ♦  Ho(1)(U 


II  il/.ll  .II  ..II  .  ,ll 

r|M  A^  cos  X  a  ♦  A^  sin  X  i 


In  Axlel  Direction  In  the  Reinforcement: 

^(r,g)  *  Jo(^r)(,[  C0*  V1*  ♦  *2  ,ln  M1*]  *  ♦  |JD(^1  + 

U 

b3r  (B3  co#  ^I#  *  B4  ,ln  ^I* 


In  Axlel  Direction  In  the  Matrix: 


«?M  -  T[^lt') ♦  b“  ho1VM)  (•,“  -  *u-  - 


»“  .In  wni)  «  +(j0(lu“r)  ♦  b“  H<°(  4iUr)  ♦ 


II  X 

COB  |i  B  ♦ 


»“  .in  m“. 


In  Tangent  lei  Direction  In  the  Reinforcement: 

$2^>s)  "  rJcJ  +  C2  *J  +  JT  Jj|lv*r|jcJ  co*  v*8  *  C4  ,ln  yI,J 


(121) 


(122) 


(123) 


(124) 
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In  Tangential  Direction  in  the  Matrix: 


Application  of  the  Solutions  to  t h*  l^rv  Condltlona  Enuatlona  for 
the  Integration  Constanta 


Prlaarlly,  the  solution*  aust  satisfy  the  boundary  conditions  expressed 
in  equations  (60)  through  (62).  The  right-hand  side  of  equations  (60) 
and  (62)  can  be  expressed  by  a  Fourier  series  in  the  dosaln 


*«***♦*« 
o  o 


The  following  can  be  stated: 


and 


(126) 


(127) 
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so  that  •  qua t Ion  (60)  can  ba  written  in  the  following  fora: 


Also,  aquation  (62)  can  ba  written  in  similar  form: 

*1^.1*  •  *'l  •  •!>  •  >1]  •  •S'kl'  •  *"l  -I1  -  ,nl] 


iitU 


The  following  is  obtained  for  aquation  (128)  by  introducing  the 
analytical  values  for  the  dlsplacesMnts: 


•ol1  *  ®x)) ^  -»(Afr  •  ’ll* 

aJ  »lnp«(l  -  B^ft  t(l  -  B1)  +  jjju1 


3^ 


(128) 


(129) 


Xha  following  la  obtained  for  th«  boundary  condition  expreaaad  by 
aquation  (129): 

Hok  *o(l  -  e1))^  «.(w1«(i .  a1)]  ♦ 


b 


•a1  •‘“(/•I1  *  e1))}  *(i  •  i1)  ♦  {^(m1  «e(i  -  o1)  «■ 
l  \>ll  -  »x)  Ji|^  *Jl '  ^j}  (»j  “'(m1*  |i  •  91))  ♦ 


* 


(130) 
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.“„«>(*“  .,(!-,«))  +b«.0(: 


\ 


TT 


i#i  .  Bxi)  ,ia  ^ 

n-l 


Squat ions  (130)  «nd  (131)  ira  identically  satisfied  only  vh«n 


further, 


Z 

1 


II  -  0 
2 


0 


la  also  valid. 
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(131) 


(132) 


(133) 


r 


I 


The  following  can  also  be  concluded: 


x1! 

|i  -  a1) 

.  HJX. 

(134) 

x11 

[i  -  a11 

1  .  an 

(135) 

A 

i-a1) 

•  f 

O 

(136) 

ill 

U  1 

1-BU 

\  .  an 

1  lo 

(137) 

where  n  ■  1,  2,  3,  4  -  -  - 

From  this  the  following  elgenvaluce  will  exlet. 

X1  - 

,1  . 

an 

(138) 

*  a1) 

.11  . 

II  r 

(139) 

M  " 

f 

Then  froa  equation  (130), 

-  e1) 

•  *“  l 

bh  y + 

*“  Hi<i)(i,m  y + *“  *0ii  -  e11)  ■>o(inT'  y + 
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(140) 


tha  first  equation  for  the  integration  constants.  Ihs  sscond  aquation 
for  tha  integration  constants  follows  from  aquation  (131): 


By  considering  tha  boundary  condition  expressed  by  equation  (61), 

•.I1  - sI)  (cf  ♦  cHi  -  »xl) + 

Jl(tvI  *.(l  -  S1))  (e3  ‘I1  -  B1))  + 


cj  .In  v1  »(l 

•  -  j*o(l 

■n 

r11  1  1 

|(C“  +  «?  ‘I1  • 

Al\\ 

1 

9  )) 

r 
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When  equation  (142)  Is  considered,  then,  all  coefficients  from  sine 
and  cosine  functions  froa  s  disappear  Identically 


0 


c“  -  c“  -  c“  -  e“  ’  0 


Equations  (120  through  (125)  can  now  be  written,  representing  the 
distortion  as  follows: 


Displacement  in  the  Reinforcement  (Radial  Direction): 


* 
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0  (142) 

(143) 

(144) 


(145) 


yi 


Displacement  in  the  Matrix  v1ladl«l  Direction): 


rll 


•r,‘)  *  ^A31{Ji(17JVTt 


+  -11  u(l)  1  — M _ JL  ]  + 

2  1  l  - 111  <J 


j4‘  t) 


ii 
*4  r 


H(l)fl  — “rr  f-fl  co.  ■  nFl  ■ 

0  l  1-1  IIloU  til  -  »“] 


Displacement  In  the  Reinforcement  (Tangential  Direction): 


«5M 


Dlsplaceawnt  In  tl.e  Matrix  (Tangential  Direction): 


*?M  ■  0 


Displacement  in  the  Reinforcement  (Axial  Direction) 


(146) 


(147) 


(148) 


(149) 


Displicement  In  the  Matrix  (Axial  Direction): 


h11  H<1> 
b2  Ho 


1  -  BU  *0 


.  .II  , 
H3rJl 


ML, 


1  -  D11  V 


r  H 


(1) 

1 


1 


(150) 


The  boundary  conditions  expressed  by  equations  (39)  and  (42)  ssist  now 
is  satisfied.  The  dlsplacesMnt  given  In  equations  (145)  and  (149) 
arc  substituted  Into  the  boundary  condition  expressed  In  equation  (39): 


1 

1  -  2V1 


S,;  t 
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Us inf  dlffsrsntlatlon  formal* •  of  ths  cyllndsr  functions, 


In  equation  (152),  all  coaffldsnts  of  ths  Bssssl  functions  JQ 
sad  r  disappear  because  the  equation  mist  be  satisfied  for 
all  r  In  notarial  Z  (reinforcement);  this  results  In  the  two 
subsequent  equations: 


Equations  (153)  and  (154)  represent  the  third  and  fourth  equations 
for  the  Integration  constants. 
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By  utilising  a  similar  approach  and  the  conclusion  that  the  expressions 
ssist  vanish  for  all  r  In  suterlal  11  (matrix),  the  subsequent  equations 
are  obtained  for  the  Integration  constants  by  Introducing  equations  (146) 
and  (ISO)  into  the  boundary  condition  expressed  in  equation  (42): 


-  0 


(155) 


i*  •  *“i  *;■  • 


(156) 


a11  A11 
%  *3 


0  (157) 


(l  *  v11)  b“  b“  -  vXI  i  AXI  -  0 


(15B) 


lie  end  conditions  of  a  finite  fiber  resulted  in  equations  (153) 
through  (158).  Mote  that  the  boundary  conditions  expressed  by 
equations  (41)  and  (44)  are  autoemtlcally  satisfied  and  hence  do  not 
furnish  equations  for  integration  constants. 

Two  additional  equations  for  the  integration  constants  are  obtained 
when  the  boundary  condition  equations  (22)  to  (24),  which  states 
that  the  stress  vanishes  on  t.. .  surface  of  the  resin,  are  satisfied. 
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Introduction  of  the  corresponding  solutions  for  distortion  from 
equations  (146)  and  (ISO)  Into  equation  (22)  yields  a  mathematical 
expression  containing  terms,  with  a  common  coefficient 


cos  nrr 


fines  the  expression  suet  be  sero  for  all  values  of  s  ,  and  since  the 
cosine  Is  not  of  this  nature,  the  complete  expression  must  vanish. 

Use  of  this  conclusion  results  In  a  ninth  equation  for  the  Integration 
constants 
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0  (159) 


A  similar  aquation  Is  obtained,  utilising  similar  conclusions,  by  using 
ths  boundary  condition  expressed  in  aquation  (14)  and  tha  solutions  from 
equations  (14*)  and  (150).  As  In  ths  former  casts,  ths  coefficients  of 
ths  trlgonoamtrlc  functions  must  disappear,  and  the  following  equation 
Is  obtained: 


.«  uO> 

*2  M1 


Two  additional  equations  for  ths  coefficient  are  obtained  from  the 
boundary  conditions  given  In  equations  (79)  and  (81).  Substituting 
oj,  and  Into  equation  (79)  yields  an  equation  containing 

cylindrical  functions  multiplied  by  cos  mr  ^  .  Since  the  equation 
base  must  be  satisfied  for  any 

cos  nrr  “ 

m 
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the  coefficient*  of  this  express  ,  jn  mu  vanish,  end  the  eleventh 
equation  for  the  coefficients  is  obtained.  Equations  (145),  (146), 
(149),  and  (150  are  utilised. 


Using  Che  boundary  condidon  expressed  In  equation  (81),  the  following 
twelfth  equation  for  the  integration  constants  is  obtained  after  intro* 
ducing  the  functions  for  o^  and  a**  : 
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The  twelve  equations  *  equations  (140),  (141),  (153),  (154),  and  (155) 
through  (162)  -  exhaust  all  technically  meaningful  boundary  conditions, 
The  coefficients  to  be  determined  by  these  equations  are 


K1  •  A1  a1 
*3  •  a3  *3 


II 


II  -II 


,11  _II 


,11  .11 


;  Aj  i  Aj  »  Aj  »3  •  A3 


' »“  .  »ll  «? 


The  twelve  equations  are  lnhosngeneous  and  linear. 

At  this  noint.  the  Physical  nrohl—  of  <ma  f lnlts-lcnsth  fiber  In  a 
matrix  la  solve*. 

The  determination  of  the  constant  Is  a  problem  of  pure  smthematlcal 
analysis.  Further  Interrelation  to  the  physical  problem  Is  unnecessary. 


Bitifliutlfla  of  tha  Intaeratlon  Constants 

because  of  equations  (153)  through  (158),  the  coefficients 


»X  b1 
*4  °3 


.II  .II  .II  .II  .II 
*4  •  B4  b3  *  *4  b2 


and 


can  be  expressed  by  the  corresponding  A  coefficients.  Equations  (141)., 
(159),  (160),  (161),  and  (162)  can  be  transformed  by  this  relation. 
Equation  (140)  can  be  taken  without  transformation  and  restated  as 


•  2  •„(#*  •  BXI)(-  l)"*1  (163) 

taly  the  results  of  an  extensive  analysis  are  given  In  the  following 
five  equations. 

Equation  (141)  becoews 


2  lo  f-  n“*1 

tt  n 


064) 
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Another  equation  It  obtained  through  modification  of  equation  (159)  by 
Introducing  equation!  (1S5)  through  (158) 


Next,  equation  (160)  la  tranaforroed  by  Introducing  equatlona  (155)  to 
(158).  Conaequently, 


Equation  (161)  beconea,  with  equatlona  (141),  (153),  and  (154) 
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finally,  aquation  (16?.)  nuat  be  traits  formed  by  expressing  all  B's  by 
A's  by  means  of  equations  (1S3)  through  (138). 
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To  s iaplify,  the  six  constant*  resulting  froa  equations  (163)  through 
(168)  can  be  defined  by 


52 


The  coefficients  C A  *o  ere  reel  nuabers.  Consequently,  the  solu¬ 
tions  expressed  in  equations  (145)  through  (150)  cen  be  written  by 
eddltlonslly  teklng  into  eccount  for  end  5^*  equetlone  (153) 
through  (158)  in  the  following  fora: 

Kelnforceaent  Distortion  in  the  Radial  Direction: 


Metrlx  Distortion  in  the  Kediel  Direction: 


telnforceaent  Distortion  in  the  Tangential  Direction: 

•1M  *  0 


(172) 
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Matrix  Distortion  in  ths  Tsngsntlsl  Direction: 


-  0 


Reinforcement  Distortion  in  ths  Axlsl  Direction: 


1  - 


lol 


i  J.li 


sin 


Matrix  Distortion  in  the  Axial  Direction: 


c*  +  e« 


r  1 
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(173) 


(174) 


(175) 


"  ? 


The  Bessel  and  Henkel  function* 


J<,M  - 1  hM 


1  H*l*|ix)  and  H*l*jlx) 


■uat  b«  recognised  ae  real  function!  for  x  >  0  .  These  functions  are 
tabulated  in  the  literature.  The  following  Identities  interrelate  the 
above  functions  and  the  so-called  Modified  Bessel  functions:  (Ref.  3) 


Especially 


•  ^'aM  ■  hi*) 

+  flOl 

2  «a1)(‘*|  *  KI-)  <l76> 

joM  •  M*) 

.Hi 

•  2  Jllu)  •  * 1  •'ll1*)  ■  hi*) 

1  "ol)M  *  SM-) 

‘2  Mia)( l*)  ■  -••ia)(1*)  •  nhl*)  <l”> 


Introducing  further  abbreviations, 


ott  a 


(178) 


S3 


(179) 


ah  b 


*2 


(l  +  4  -  2vXJ  *l 

jl  ♦  vn)(l  -  2V11)  ‘  "* 

In  Addition,  th«  following  is  sst: 


(180) 


(181) 


Ths  following  equation  exists  for  the  determination  of  the 
C4  (i  •  1  to  6) 


where 


(1  to  6) 


The  36  B's  are  defined  according  to  the  equations  used  to  this  point  of 
development: 


(182) 


(183) 


B13  "  kj’  Xl(kl) 
B14  *  ^*l(kl) 


(184) 

(185) 


(136) 

(187) 

(188) 

(189) 

(190) 

(191) 

(192) 

(193) 
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(194) 

(195) 

(196) 

(197) 

(198) 

(195) 

(200) 

(201) 

(202) 

(203) 

(204) 

(205) 
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(206) 

(207) 

(208) 

(209) 

(210) 

(211) 

(212) 

(213) 

(214) 

(215) 
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Hum r leal  Determination  of  C’s: 
Gl*en  dimensions: 


(216) 


(217) 


(218) 

(219) 

(220) 

(221) 

(222) 

(223) 
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Basic  computations: 


(224) 

(225) 

(226) 

(227) 


(228) 


T 


I 


(255) 

(256) 

(257) 

(258) 

(259) 

(260) 

(261) 

(262) 

(263) 

(264) 
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(265) 


(266) 


(267) 


(268) 


(269) 


(270) 


m  run  a  A  MM  rYLimn 

The  basic  equations  (1)  through  (10)  are,  In  this  case,  also  valid,  while 
the  boundary  conditions  change  Insofar  as  there  Is  no  finite  length 
expression  given.  After  separation  of  the  variables  by  Introducing  a 
Bernoulli  product  of  the  fora 

5  •  *<r)  «tU* 

and  further  Introducing  In  the  obtained  differential  equation  the 
potential 


$(*) 


.2  x  dx 
dx 


K 
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where  x 
obtained 


Xr  and  X  It  the  eigenvalue.  The  following  solutions  art 


t[(r,sj  -  A^  4  JojlXrJ  ♦  *2i  Jl(i^r)j  ‘  co#  *273* 

*“|r,sj  -  A2r  4  ^4  Jljir  H<l)jiXr)  4  »4  H<l)|iXr)J  cos  Xs  (274) 

(2|r,s j  -  A4  4  J»5  JojiXr|  4  B6ir  J^iXrjj  cos  Xs  (275) 

«“|r,s)  •  A5  4  f g?i  H^l)(  iXr  j  4  §8r  hJL>J  iXr|l  coa  Xs  (276) 


Applying  tha  solutions  to  ths  corresponding  boundary  conditions ,  sight 
equations  for  ths  Integration  constants  B's  result.  His  system  deter¬ 
minant  was  homogeneous  and  the  following  four  eigenvalue  equations  were 
obtained; 


J*|xb)  +  Jl(Xb)  "  0  (277) 

Hll)  (Xb)  +  Hol)  (Xb)  *  0  (278) 

+  -i  J.M  JiM  • 0  U7,) 

’o1)2|H  +  Ml0)Vb)  '  »  Hol)(Xb)  "l(l)M  •  0  (MW 
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tfr 

The  SMthods  given  In  Computation  of  Kankel  Functions  did  not  result  in 
conjugated  coaplex  roots.  This  indicates  that  all  B's  are  sero  in  the 
solutions  expressed  in  equations  (273)  through  (276),  and  the  stresses 
for  the  Infinite  fiber  becosM 


*  Coaoutetlon  of  Henkel  Functions.  National  Bureau  of  Standards  Report 
No.  216. 
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•  2vI°n(H 


(286) 


,  rII  V11 

_ 2E _ __ 

V  -  a11  ♦ 

P 

M 

1 

P 

M 

K 

-  - 

MM  l  +  iijx-av1® 

L 

♦  V1 

E 

1  +  v11 

K) 

(287) 


a 


I 

12 


II 

°13 


(288) 


Since  Che  infinite  undistrubed  fiber  is  not  realizable,  equations  (280) 
through  (288)  are  values  which  ars  difficult  to  verify  by  test;  however, 
they  can  be  used  to  provide  an  estimate  of  the  magnitude  of  residual 
stresses. 

For  Instance,  if  the  relation  between  the  reinforcement  and  resin 
modulus  were 


•  v2  ■  1/4  and  b  »  a 

and  assuming  that  the  total  difference  of  contraction  Including  polym¬ 
erisation  were  4%,  then  by  using  equation  (281)  the  following  would 
exist: 


This  indicates  that  the  radial  compression  of  a  glass  fiber  with  a 

6 

modulus  of  11  x  10  psl  would  be  15,700  psl  on  the  surface  of  the  fiber 
(r  ■  a)  . 
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For  b  ■  10  *  ,  the  axial,  tensile  stress  In  the  resin  would  be 

14,000  psl  which  would  cause  cracking;  the  resin  would  therefore  become 
a  finite  length  and  the  equation  for  Infinite  length  would  not  be 
adequate.  For  this  reason,  no  more  emphasis  Is  given  to  the  Infinite 
length  fiber  with  assumed  boundary  conditions. 

The  problem  would  be  different  If  the  boundary  conditions  that  are  no* 
phvslcally  Imposed  were  derived  In  the  fora  of  a  differential  equation 
which  provides  for  cracking  In  certain  unknown  distances. 

ihe  Magma  mam 

Assume  that  an  Infinite  number  of  fibers  are  packed  so  that  In  the 
limiting  case,  the  tightest  packing  Is  possible.  A  configuration  such 
as  this  would  possess  a  certain  synsaetry  which  would  allow  for  rigorous 
analytical  treatment. 

Uwlng  the  configuration  depicted  In  Figure  3,  any  one  fiber  can  be  con¬ 
sidered  as  the  central  fiber  surrounded  by  six  others  (all  fibers  are 
Imbedded  In  the  matrix) . 

Boundary  C  >ndltloni  for  the  Multifiber  Problem 

By  assuming  three-dimensional  shrinkage,  a  boundary  condition  Is  obtained 
on  the  Interface  between  fiber  and  resin  which  Is  similar  to  the  boundary 
condition  for  a  single  fiber  In  a  resin  cylinder,  the  only  difference 
being  that  the  distortions  will  depend  also  on  the  angle  OP 

n  [sl1  ’  eJ»  9  *  *(l  ■  9l)j  "  5“[*o(l  "  9  •  *ll  ‘  )] 

•  *0\ex  -  32J  (289) 
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0  (250) 


*2  |*o  I1  "  8l)  •  9  ’  ‘l1  ‘  9ll]  '  ^[‘ol1  *  ’  9  1  *1'  '  ^2)  * 

«3  [*oll  -  ®i)  •  9  •  ‘I1  •  »i)]  -  ‘“[‘cl1  -  e2|  •  9  •  ‘I1  •  82)] 

•  «(b1  -  »,)  (291) 

Where  5^  it  the  displacement  in  direction  1  (or  r)  in  the  reinforce* 
ment,  Is  the  radius  of  the  fiber  before  distortion  through  shrinkages 
takes  place,  and  z  is,  similarly,  the  coordinate  before  distortion. 

Referring  to  Figure  3,  it  can  be  assusted  for  reasons  of  syaesetry  that 
the  hcxagonals  will  remain  regular  during  the  shrinking  process,  and  will 
become  proportionately  smaller  as  the  total  composite  shrinks. 

The  displacement  vector  perpendicular  to  OB  must  vanish  if  the  triangle 
sites  OAB  remain  straight. 

For  the  reinforcement: 

[r(l  -  a1)  .  0  ,  >{l  -  6*jJ  -  0 

For  the  matrix: 

[r|l  -  e1)  .  0  .  <|l  -  s11)!  -  0 

For  the  OA  line: 

In  the  reinforcement 

*2  [rll  *  0l)  •  6  *  *ll  ’  6l)]  '  0 


(252) 


(253) 


(294) 
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and  in  the  Matrix: 


['  (>  -  f>“|  ’  0 

Along  tha  AB  line: 


(295) 


M 


and 


Kquatlona  (292)  through  (295)  nuat  be  valid  for  all  valuaa  of  r  and  z  , 
ana  aquation  (296)  for  all  valuaa  of  $  In  tha  region  where 

0  a  qp  a  J 

Since  the  triangle  QAB  repeats  12  tinea  in  each  hexagonal,  It  is  possible 
to  solve  the  problcn  by  considering  QAB  only. 
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Differential  Equation!  for  Dlsplacea inf 

General  partial  differential  equations  for  distortions  considering  the 

* 

cylindrical  coordinate  system  were  derived.  These  equations 


a2?,  .  »2«.  i  a2*,  i  a2t> 


i  a2{,  ,  a2t, 

1  -  +  2  r2  a<*>2  +  2  9.2  +  2<2  -  2v>  '  Sr3* 


a2« 


_  3  *  4v 


2(1  -  2v)  Sras  2(1  -  2v)  f2  a$ 


Lix  1 

1  -  2v 


i 

r  ar  r 


0  (297) 


a2«, 


1  a  s!h+i  !!*2+. 

2  Sr2  1  •  r2  9<p2  2  9.2  2<l  '  2v>  r  8r» 


Jl 


1 


a2!, 


2(1  -  2v)  r  a<p  a* 


I  1  f  3  -  4v  i  +  .  ii 

2r  1  -  2v  r  dtp  ar  r 


•  0  (298) 


1  ?!*1  +  _L. 
2.2,2 
ar  2  r 


i!b+ 

s»2  1  - 2v  9.2 


JL 


a2l, 


2(1  -  2v)  ara* 


2(1  -  2v)  r  depds  +  2(1  -  2v)  r  a*  *  2r  ar 


-  0 


(299) 


*  Equation  (A73)  In  the  Appendix. 
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The  following  identities  ere  noted. 


X 


2(1  -  2v) 


1  -  v  .  1 

1  -  2v  2 


JLZ-toL  ,  L JLJd  +  1 

2(1  -  2v)  1  -  2v  2 


These  ere  used  In  equetions  (297)  through  (299),  which  on  reerrengement 
end  slmpliflcetion  become  the  following  sequential  equetions: 


L.:,  3L 

x[i 

£!■ «,)  +£  + 

i  2 

1  -  2v 

dr  [r 

r  2 

X 

X  [r 

|2*i.2!i]l  +  i 

xr 

2r 

ds  [ 

dr  ]J  2 

d!p[ 

l  -  V 

1  X 

l  2*i  +  2*i  +  i 

XI 

1  -  2v 

r  d5P 

r  d$  d*  r 

dr  \ 

l  2*i 

.2  dCP 


2  ^2 


1 

r  dr 


1- 


0  (300) 


1  X 

2  dr 


(r  [*  Ir  *  ^D+  2  ^[Ir1  *  f  5^] 


0  (301) 


j^_*.  .a-fi  -a.ir  t » +2ii  +  i  2!il 

1  -  2v  d*  [r  dr  |r  *1/  d*  r  dQP  J 


X 

X 

+  xxfi 

2fal 

2r 

dr 

r\dr 

ds 

,  2r  dQP  Lr 

dCP 

3.  J 

-  0  (302) 


Some  recurring  functions  in  the  previous  equetions  ere  noted  end  defined 
es  F's  es  follows. 


LUL  1  X 

1  -  2v  i 


£lr  {i| 


(303) 
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(304) 


L 

2 


2  [a*  r  dcp  '4 


(305) 


(306) 


Using  equations  (303)  through  (306)  in  equations  (300)  through  (302) 
yields  the  following  partial  differential  equations: 


sf.  ,  ar  ar, 

—i.  +  -  -—4  4  —i  .  o 

ar  r  a*  dCP 


(307) 


1 

r 


0 


(308) 


r  a  Cp 


(309) 


Now  F's  can  be  considered  as  potential  functions.  First,  efforts  would 
be  sttdc  to  find  them.  In  addition  to  the  previous  three  equations,  the 
procedure  that  follows  renders  one  more  equation,  interrelating  the  F's. 

Multiplying  equation  (305)  by  r^  and  equation  (306)  by  r  yields 


2  lacp  dr 


r2  F, 


(310) 
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[illiil. !*i]  .  rF 

St  S9  J  r  r4 


(311) 


Vlartial  differentiation  of  equation  (310)  with  respect  to  a  and 
equation  (311)  with  reepecc  to  r  glues  the  following: 


2  I  dsddp  Sadr 

4fthi.il 

2  ^  5r5*  5r59J 


(312) 


i  lr  r*l 


(313) 


tlon  of  equations  (312)  and  (313)  gives 


i[i!L .  !!*i]  .  r2  ^i  +  -a-ir  r ) 

2  |5*59  drdqpl  5a  5r  |r  r4| 

Partial  differentiation  of  equation  (304)  with  respect  to  9  gives 


(314) 


l»Si  | 


.  ih] . 

I  595a  “  595rJ 
I  5*59  ’  5r59l 


1 

r  59 


(315) 


Noting  that  left-hand  sides  of  equations  (314)  and  (315)  are  identical, 


one  obtains 


A  ^  2  5^  5(r  F4 

r  59  r  5*  ’  5r 


(316) 
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Separation  of  the  PotintUU: 

V 

Multiplying  equation  (308)  by  r  and  partially  differentiating  with 
reapect  to  s  ,  one  obtains 


r 


dsdr  d«d9 


(317) 


fbrtial  differentiation  of  equation  (309)  with  reapect  to  cp  gives 


1 

r 


-  0 


SuMation  of  equations  (317)  and  (318)  renders 


i  •  e  • , 


t  \ 

8s2  *  r  a*2  ’  &*3r  '  apdslr  'll  "  ° 


Froai  equation  (307), 


fcP 


8F.  .  8F 

rrk  l  1  r- 

dr  r  ds 


(318) 


(319) 


(320) 
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Substitution  of 


dr 

dtp 


In  equation  (319)  renders 


r 


±1,1  ±i.k L  +  ±ft  + 

».1  r  w2  8,dr  9r  L  9r  J*  J 


Division  by  r  gives 


(321) 


fkrtlal  differentiation  of  equation  (307)  with  respect  to  s  gives 


d2FI  ^  dV^  d2F3 
dtdr  +  r  Sf2  *  d*d«P 


(322) 


Flsrtlal  differentiation  of  equation  (309)  with  respect  to  r  gives 


Subtraction  of  equation  (323)  fron  equation  (322)  renders 


1 

l  !IiL  ^  +  .a. 

s 

1 

5!il  . 

r 

».2  9r  l» 

*  dr  J  dsdcp  dr 

r 

dtp  J 

(323) 
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l.e 


» 


Frotn  equation  (316),  rearrengeaent  give* 

5!i  .  x  !Iz .  x  ±jA 

da  f3  dCP  f2  dr 


(324) 


(325) 


&F. 

Substituting  for  froa  aquation  (325)  In  equation  (324)  and 

multiplying  the  resulting  expression  by  r  gives 


!^+r  xli  !!i)  +  x 

J.2  Jr  lr  Jr  J  ,2  ^2 


81 


T” 


026) 


-  -1 

ft.  d  F. 

-i  +  — i  - 

d<P 

L  r  dr 

dr  r  r4j 

r  d$> 


From  equation  (309), 


r  Sep  3*  r  Jr 


The  equation  (326),  after  rearrangement,  become* 


Mow 


a  i  fiLi  Hi 

dr  |  r  dr  /  r  dr 


(327) 
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o 


dF 

— -»  -f 

d*F2 

dr 

dr 

±[r 

dr2 

dr  lr 

dr 

.  1 


&r. 


r  I  dr 


+  r 


dY 


dr 


(328) 


Use  of  equation  (328)  in  equation  (328)  gives 


d2r, 


,  d2P. 


-  2 


d* 


(329) 


Multiplying  equation  (308)  by  r  and  then  partially  differentiating  the 
reeulting  expression  with  respect  to  r  gives 


drd5P 


.  ± 

dr 


fr  ^  |r  F-] 

[  dr  V  3| 

J 

(330) 


Bsrtial  differentiation  of  equation  (316)  with  respect  to  s  renders 


r  d*dcp 


-  t 


(331) 


Sunmation  of  equations  (330)  and  (331)  gives 


dcp 


(332) 
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From  equation  (307), 


ar  ^  ar^  ar. 

ar  +  r  a*  1  a«p 


Then  equation  (332),  after  multiply in*  by 


*nd  rearranging,  become* 
r 


aV, 


a<p 


X 

2 

r 


*[' 


« 

FjJ  0 


(333) 


Partial  differentiation  of  equation  (308)  with  reapect  to  r.  glvea 


aadQPl  r 


-  0 


(334) 


Dividing  equation  (309)  by  -  r  and  partially  differentiating  the 
reaultlng  expreaalon  with  reapect  to  9  ,  one  obtalna 


jl[x 

avae  i  r  /  acp  l  r2 


(335) 


Summation  of  equatlona  (334)  and  (335)  rendera 


From  equation  (316), 


1 

r 


Multiplying  by  r  and  martially  differentiating  with  respect  to  x 


j:r 


arfcflp 


Dividing  by  r  gives 


£['3  +  ^  [r  £lrf*l] 


aFa  ilil  .  a  ±[tJ  +  Ji.[r  -a./rF|l 
dQplr2  hT  J  r2  dr  L  3*  J  r2  &r  L  dr  * 


Substituting  for 


a[a 

ftl 

»L2 

9r  J 

from  equation  (337) 

In  equation  (336) 

gives 

A  A.  [r3  fll  . 

a2  ,  «  i2 

>2 

F  ♦  *— 

(!i]  + 

r2  ir  Lr  »•  J 

dr da  1  3/  jf2 

4  a*2 

lr2J 

i-JLf 

r  -4- 

2  dr  1 

df 

r  l 
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(337) 


or,  rearranging 


f2  dz  a* 


Separation  of  the  Varlablea 

Separation  ol  {'a  in  the  Form  of  Partial  Differential  Equation 
Containing  F'a: 

Reference  la  made  to  equationa  (303)  through  (306). 

Separation  of  ,  the  dlatortlon  In  radial  directlona: 

From  equation  (304), 


(338) 


(339) 
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From  equation  (305), 


A. 

8r 


2 

r 


Multiplying  equation  (303)  by  ^  gives 


-*■  fr  *  ] 

dr  r  dr  \T  *l| 


i  t  8ii 

lr  5i)  +  51 

h\ 

1  - 

p 

fcp 

r  1 

|  1  •  V 

1 

respect  to  t 

&2 

fii) 

|  .  L-:  1y 

8F 

’i)J  a*ar 

acpdr 

r 

1  -  V 

&r 

Now 


(340) 


(341) 


(342) 


Using  equations  (339)  and  (340),  equation  (342)  would  becone 


?.tdr  JPPdr 


h) 

.  A 

2F/ 

.x  !*a  + 

1 

dz 

.a* 

r 

f2  &P 
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(343) 


a2«, 


a* 


+  i  if. 

2  «l1 

2  F  r  ♦  — 

r2  aU 

3  a«p  J 

1  Hi 

2 

-  -a-  — %  -  2 

ar3  .  j2^ 

r  ds 

2  acp 

r 

»  *7  w2 

Subtracting  aquation  (343)  from  equation  (341) 


.  1-1 
1  -  v 


*1 
r  a* 


ar_  ,  h 

■  ■*  +  ■*  — ^ 

acp  r2  dep 


(344) 


Separation  of  *2  »  th«  distortion  in  tangential  direction  yields: 
Prom  equation  (305), 


Prom  equation  (306), 


*  £('  +2  r'r3 

'  r  £  lr  {2)  +  2  r3  F3 


(345) 


(346) 
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Multiplying  equation  (303)  by  j—"1  ^  and  then  partially  differentiating 
with  respect  to  op  , 


L-I  2v 
1  -  v 


Multiplying  (r)  results  in 


dtp 


♦  r 


\  I  M2 


L  :.1Y 

1  -  V 


(347) 


Using  equations  (345)  and  (346)  in  equation  (347)  gives 


Rearrangement  renders 
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Separation  of  5j  ,  the  distortion  in  axial  direction  results: 
From  equation  (304), 


Fro*  equation  (306), 


dC9  ,  d( 

— 4  -  2  F  +  *  — 1 

ds  4  r  d<P 


From  equation  (303,  by  nultlplylng  Ly 


1  •  2v 

J  ^  “  r  one  obtains 


dr 


1  -  2v 
1  -  v 


Partial  differentiation  with  respect  to  z  gives 


Use  of  equations  (349)  and  (350)  in  equation  (351), 


i  ± 

r  »r 


r2  a*2 


l-v  *r 


(349) 


(350) 


(351) 


(352) 
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Rearranging 


+  1  f  r 

r  dr 


L 

.  l| 

dF2  dF\\ 

1  -  V 

a* 

r  | 

dr  +  dcp  | 

sournow  FOR  S2 

Assuming  that  the  F's  are  zeros,  equation  (348)  becomes 


(353) 


Let 


S2  -  R2(r)  02(cp)  Z2(z)  (354) 

Then  equation  (353)  becomes 

*  2  "  2  " 

R2  ^2  Z2  ♦  3r  R2  02  Z2  +  r  R2  ^2  Z2  +  r  *2  ®2  Z2  * 

R2  *2  Z2  "  0  (353) 
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Dividing  equation  (355)  by  ?2 


H  II 


^2  2  *2  2  Z2  ^2 
1  +  3r  ^  +  r  V  +  r  V  +  .  " 

*2  *2  Z2  *2 


or 


M 


r2  ^  3r  ^  +  r2  +  1 
*2  R2  Z2 


h 

*■> 


-  x: 


(356) 


where  X2  is  s  constant  independent  of  r  ,  t  ,  and  cp  . 


02  ♦  X2  02 


(357) 


and 


r*^+3r^+rJ^+l 
*2  *2  2 


♦  X 


(35S) 


The  last  equation  can  be  rewritten  as 


i!l  +  2  I2  +  .L  i 

I,  r  I  2  1 

2  r 


-  X 


h 

h 


(359) 


where  n  is  a  general  constant  Independent  of  both  r  and  z 
equation  (359)  nay  be  separated  as  follows: 


"  2 
Z2  4  »  ZJ 


(360) 
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I 


and 


•I 


+  i»' 

r  2 


(361) 


The  general  eolation  of  equation  (357)  la 

$2 (CP)  •  coa  |X2*>  ♦  cpQ|  (362) 


where  qpQ  la  the  phaae  angle  and  B^  an  arbitrary  conatant. 

In  the  caae  of  alx  ayueetrlcally  apaced  flbera  aurrounding  a  central  fiber 
it  la  expected  that 


02(cp) 


0, 


where  J  la  any  Integer. 
Thla  »ean a  that 


coa  \,QP  ♦  cp 
2  o 


For  equation  (364)  to  be  aatlafled  for  any  m  , 


2  tt  k 


(363) 


(364) 
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(365) 


i.e. , 


X2  -  6  k 


vhara  k  la  Any  Intend. 


Tharafora,  aquation  (362)  becoaii 


02  -  cos  1 6  ktp  +  5PoJ 


(366) 


and  aquation  (361)  bacoaas,  aftar  dividing  by  n  , 


i  *%  a 

m!  ir  ~2l  uv 


(367) 


Lat 


lit 


(368) 


Than 


<UL. 
~2  dr 


dft 

__  za 

2  dx 


4ft 

dr 


dl 

Zi 

dx 


(iu) 


d». 

lrjj  dx 


.  i  fa 

X  dx 


(369) 


± 

2.2 
li  dr 


.  ±  -i-ffi  4ft  \  4ft 

2  dx  Idx  dr 
u  >• 


l\  4 ft 

rj  dt 


.  (li^2 

M2  dx2 


dV 


dx 


(370) 
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Using  equations  (368),  (369),  and  (370)  In  equaelon  (367),  and  Multiplying 
the  resulting  expression  by  -  1  ,  yields 


dx 


(371) 


To  transfona  the  last  equation  to  a  Bessel  fora. 


t2  ■  x°  y  (x) 


Then 


(372) 


dx 


0  x®”1  y(x)  ♦  *fl  y'(x) 


(373) 


dV 


dx 


0  (0-1)  x9"2  v(x)  ♦  20  x9*1  y'(x)  ♦  xa  y"(x) 


(374) 


Using  equations  (372)  through  (374)  In  equation  (371)  and  dividing  the 

a 

resulting  expression  by  x  results  in 

♦  y'  (x)  ♦  |i  -  ^.,1  2*  ijy  ♦  {3a  ♦  a  (a-D }  " 


y"  (x) 


Rewriting  equation  (373), 


0  (375) 


y"(x) 


+  1Z1-1  y.(x)  +  £  .  36  k?  •  i  -  136  +  fl 


1  (376) 
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6  la  to  ba  choffn  so  that 


i(C<  i 


20+3  -  1 


6  -  -  1 


Th«n  aquation  (376)  bacoaaa 


x 


(377) 


Equation  (377)  la  a  Bataal  aquation  and 


y  -  »2  J6k(x)  ♦  »,  SV<x) 


(378) 


Ihls  la  tha  gancral  aolutlon  vhara  and  ara  arbitrary  conatanta. 


from  aquation  (372). 


■  _ 

*2  *  y 


Slnca  0  ■  -  l  , 


R,  -  * 

2  x 


Subatltutlnt  for  y  froa  aquation  (370) 


(379) 
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Became  x  •  Ipr  ,  equation  (379)  la 


*2<r) 


pr  L  1  6k 


§ 

(lyr)  +  -j1  (lyr) 

i 


The  general  aolutlon  of  equation  (360)  la 


*2  *  »4  •!«  |  M«  ♦  *0) 


where  and  *Q  are  arbitrary  conatanta. 

Subatltutlng  for  *2(r)  »  02^  »  and  *2(«)  fro«  equatlona  (380), 
(366),  and  (381)  reapectlvely  In  equation  (334) 


(2(r,  cp,  a) 


[*1  co*  I6**  +  ®o)][*4  •lB  I  M»  ♦  *o)] 

£  [»5  +  »*  <JW>]  * 

coa  1 6kcp  +  #1  x  a  in  (  pa  ♦  a  j 


where 


B1  B4  »2 


and 


B 


€ 


fjA h 

i 


(380) 


(381) 


(382) 
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Note  that  u  may  assume  many  values,  dependent  upon  boundary  conditions. 

Let  ^  ^  where  n  varies  from  1  to  infinity  and  is  an  Integer. 

As  stated  earlier,  k  varies  from  0  to  infinity  and  is  an  integer.  A 

general  assumption  is  that  the  constants  Be  ,  B-  ,  5p  ,  and  x  and 

j  o  o  o 

corresponding  *-n  equation  (362)  would  be  different  for  different 

combinations  of  k  and  n  .  To  indicate  this,  the  said  parameters  would 
be  properly  subscripted.  Then  equation  (382)  would  take  the  following 
form. 

"ink  J6k(S,r)  +  86„k  ■tfJWl]  X 

c.  (*W+  «Ponk|  *  ♦  «onk) 


?2nk(r’ 


•  *2nk»2nkZ2nk  <383> 


Mow  the  total  distortion  of  ^  i#  * 
by  equation  (383),  or 


tion  of  all  ?2nk  given 


t,(r,  cp,  s) 


■tz 


*2nk(r'  «’•  *> 


m  m 

■  E  E  v  [B5nk  j«>Kr) + '‘a*  vi] 


n-l  k-0 


[co‘  I6*®  +  "onk )]  x  [,ln  ( V  +  *onk)  (384) 
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SOLUTION  FOR  g3 


Consider  equation  (3S2).  F's  are  to  be  taken  as  tero. 


i  •  e  •  , 


0 


Let 


«3<r)  03(V)  Z3(r) 


Then  equation  (386)  becomes,  after  dividing  by  and  rearranging 


R3  Zj  ♦  -3 


Z„  ♦  1  r! 


r  3  ^3  Z3 


2  Z3  R3  *3 


i.e 


•  # 


ft  H  •  |« 

fa  .  ia  .  i  li4j.  fa 

z3  R3  r  R3  r2  ^3 


0 


••  •  11 

fa  + 1  !a4.L  fa 

R3  r  R3  r2  ®3 


where  w  is  e  constant  eigenvalue,  Independent  of  r  ,  $  ,  or  s  . 


(385) 


(386) 


(387) 


(388) 
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Equation  (388)  on  separation  gives  the  following  two  equations. 


ii 


♦  w 


0 


and 


•• 


♦ 


i 

r 


M 


2 

w 


Equation  (390)  Is 


ii 


-a+r-i- 


2  2 
w  r 


(389) 


(390) 


(391) 


where  la  a  constant.  Equation  (391),  on  separation,  gives 


«3  +  xj  #3  •  0 


(392) 


and 


2  ^3  I  22  2 \ 


(393) 


Let 


x  -  Irw 


(394) 


Then 


dR. 

dr 


~  dR*  j  v  dR.  dR. 

*  ii  ■  i  ■  (iw)  ■  x  ^ 

lw  dx  dr  lw  dx  J  dx 


(395) 
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Similarly 


2 

r 


Then  aquation  (393)  becomes 


2 

x 


(396) 


(397) 


Equation  (397)  ia  a  Bessel  equation  and  its  general  solution  is 


*3(x) 


ci  v° +  Cj  Hs(x) 


(398) 


A  similar  argument  can  be  utilised  in  arriving  at  a  solution,  from 
equation  (366),  for  aquation  (337).  Then  one  obtains,  with  \ ^  •  6g  , 

a  solution  for  equation  (392). 

*3  *  C3  cos  (  6gqp  ♦  5pJ  (399) 

where  g  is  any  Integer  and  ?0  an  arbitrary  constant. 

Noting  that  x  -  irw  and  \ ^  •  6g  ,  equation  (398)  becosws 

R3(r)  -  Cx  J6g(irw)  +  C2  H^irw)  (400) 

The  general  solution  of  equation  (389)  is 

Z^U)  ■  sin  jwt  ♦  sqJ  (401) 

where  zQ  la  an  arbitrary  constant. 
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Substituting  for  R^(r)  ,  ^(ep)  ,  and  Z^(s)  fron  aquation*  (400), 
(399),  and  (401)  raapactivaly  In  aquation  (387) 


?3(r.AM) 


[■ 


Virv) 


+  C6  H^drwjjJeo. 


(6W  +  «io)[.ln(w«  +  tj 


(402) 


where 

c  •  c  c  c 

U1  U3  U4 


■ote  that  v  My  assua«  Many  valuas,  dapandlng  upon  tha  boundary  conditions. 
Lat  w  •  where  m  varlas  froM  1  to  infinity  and  is  an  Intagar.  Tha 
constants  Cj  ,  ,  «po  ,  and  sq  and  corraspondlng  ?3  would  ba  ,  In 

general,  dlffarant  for  dlffarant  conblnatlono  of  g  and  •  .  To  indicate 
this,  propar  subscripts  would  ba  used,  giving  equation  (402)  tha  following 
form 

*3-.  *  [C3-.Jtk(l™.)  +  e.M"»i>|1"-)][e0,(‘W  + 

*o^)][,l0lv  +  \*j]  (403) 

Tha  total  (3  Is  a  suanntlon  of  all  *s  given  by  aquation  (403): 

•  •  •  • 
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♦  c 


6«g 


I6*®  +  *«sl][,lB  1 V  +  'j] 


(404) 


GOMPLMEJfTAAY  SOLUTIONS  FOt  ( 


1 


Taking  F' ■  at  saro,  aquation  (344)  ba conag, 


4ftar  alapllflcatlon,  aquation  (403)  baconta 


Conplaaantary  function: 
Conaldar 


0 


Ut 


*  ^(r)  0l(9)  ZjU) 


(405) 


(406) 


(407) 


(408) 
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be  the  solution  of  equation  (407).  Substituting  from  equation  (408)  for 
in  aquation  (407)  and  dividing  the  result  by  g^ 


(409) 


Readjusting, 


(410) 


where  0  is  constant.  Separating  equation  (410),  the  following  two 
equations  result: 


"  2 

zl  +  r  •  0 


n  •  »i 


2 

Consider  equation  (412).  Multiplication  by  r  gives 


or 


"  9 


‘i  n 


(411) 


(412) 


(413) 


(414) 
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*nd 


5  it  • 

r  r.  r  R. 


,  2  2  2 
• 1  *  *  r  -  s  -  0 

i  i  1 


(415) 


Multiplying  tquAtlon  (415)  by 


2  ii  •  2  2  2 

r  Rj  +  r  Rx  •  RL  •  B  r  R^  -  Rj 


(416) 


Ut 


x  -  l|r 


(417) 


Th«n 


dR 

r  ^ 


dR,  . 
JL  U  di 
IB  dx  dr 


dR.  2  d^ 
2 


x  7-*  r 
dx 


Uk 


10 


2  Al«. 

2  dr  Ldr 


10 


2  .  TdR 

? 

2  dx 


al| 

x  dr  J 


dV 


2  "1  .2  2  t2  .2  2  2 

x  — ,  *  0  r  ■  1  0  r  -  x 

dx 


Then  tquAtlon  (416)  btcotttt 


(418) 
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Ths  solution  of  squstlon  (411)  Is 

Z1  •  *1  ‘ln  I  »,  +  *„)  <4,s> 

Ths  solution  of  squstlon  (414)  Is 

*1  "  A2  C0B  (  *  Vo)  (420> 

vhsrs  «o  and  cpQ  srs  arbitrary  constants. 

Rots  that 

0 1  (*)  -  j  +  ^  )  (421) 

vhsrs  J  Is  any  lntagar(  la  sxpsctsd,  dus  to  slx-flbar  syasatry  around 
ths  csntral  flbsr  Thsn  axaxlnatlon  of  squstlon  (420)  glvss 

-  ba 

vhsrs  a  !•  an  lntsgsr. 

How  aquation  (420)  bscoass 

-  A2  cos  1 6afip  ♦  5pJ 

Squstlon  (41B)  Is  a  Bssssl  aquation.  Noting  that  x  ■  10r  from 

aquation  (417)  and  ■  6k  fro*  aquation  (422),  ths  gsnsrai  solution 

of  aquation  (418)  aay  bs  rsvrlttsn  as 

I^r)  -  A3  J  - (10r)  4A4  H(l> _ <10r)  (424) 

V 36a2  +  1  V36a2  ♦  1 


(422) 


(423) 
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Combining  equation*  (408),  (419),  (423),  And  (424), 


(L(r,9,8)  -  Rt(r)  ^((p)  Z^a) 


[' 


-  r  a5  j 


,0) 


— - (18r)  +  A.  H  ' _ (igr) 

<Zj7i  vwt:  ji 


>][cot(&a9  + 


«P, 


)][“lnK +  *o)] 


(425) 


where  A.  ■  A.  A.  A.  and  A,  ■  A.  A.  A.  . 
5  312  ^  412 


th 


Let  f  be  the  p  value  of  |  .  Noting  the  dependence  of  A.  , 
P  5 

>0 

Equation  (425)  then  becoMe 


A.  ,  Cp  ,  a  and  (.  on  p  and  a  ,  "proper  eubscripte  would  be  ueed. 
boo  i 


«lPa(r-®-*)  *  %aJ<— 2 - K'l 

L  V36a  ♦  l 


A.  H 
6pa 


(l) 


_ |lBpr  |j  £co«  |6acp  + 


'opJlMV  +  ‘opal] 


Suanetion  of  all  5^^  given  by  equation  (426)  would  give  the  total 


(426) 


I 


1  ' 


^(r.qp.a)  - 


^lpa(r,v,,>  " 
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Au.  H(l) 

^  V 


l—i — llv)l[co,l 

36a  +  1  J1 


Islnfp  i  4-  g  M 

»  p  opa»J 


ropaj||  »p  opai 


Por  caapleteneee,  the  other  two  {'•  will  be  stated  herein. 


£  [»5nk  J6kM  + 

Snk  <l'M]  [~K *  *j][Hv  +  *oJ 

«3  • 

■-1  f-0 

cfa*  "tiM1”-)]  ["•t6-* + *o^)][,ibIv  +  *.».)] 

Particular  Integral  for  Jilil 

The  particular  Integral  of  equation  (406)  auat  be  found.  The  right- 
hand  side  la 
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(427) 


From  aquation  (384), 


l2<r.«M> 

n"l 


•  • 

■IS- 

n"l  k*0 


2nk  *2nk  Z2nk 


whara  J 


£[■ 


’2nk  u_r  |"5nk  "ftkl^n*)  "6nk 


♦  »< 


•CM] 


»2nk  *  “*  H  +  %k) 


(428) 


*2nk 


-  .InJ^ 


1  +  conk] 


Than  tha  right-hand  alda  of  aquation  (406)  would  ba 


■lii 

n"l  k-0 


2  *2nk  Z2nk  *2nk 


Ut 


Y,  X! - T**  ***  ,lnitWP  +  ®onk) 


n«l  k-0 


12  k  R,  . 

v  _ _ 2b)Sl 

2nk  _2 


(425) 
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I 


(430) 


and 

32nk  ’  H‘*P  +  »onk) 

Then  the  right-hand  aide  of  equation  (429)  becoaea 


■tt 

n"l  k*0 


2nk  2nk  w2nk 


Ut 

l2nk  Z2»k  *2»k  -  (rl*ht-h»nd  .l<U)ok 
Then  the  right-hand  aide  of  aquation  (432)  ia 


(431) 


(432) 


(433) 


(right-hand  aide)  . 

nk 


(434) 


To  obtain  the  deaired  particular  integral,  each  (right-hand  aide)  . 

nk 

would  be  taken  aa  right-hand  for  equation  (406).  Particular  integrala 
for  auch  equationa  would  be  found  and  all  auch  reaulting  expreaaiona 
auaned.  For  exaag>le, 


itiliak  +  -L 

ir  2 


no 


(right-hand  tide) 


(435) 


nk 


where  th«  particular  integral  for  the  ipeclflc  n,k 

combination. 


Let 


*1  nk  “  Rnk*r)  *nk(qp)  Znk(,) 


(436) 


Substituting  from  equation  (436)  in  equation  (435)  and  rearranging 
results  in 


R„k(r) 


1  0nk(5p) 

«ferR»k(r) 


*2nk  Z2»k  »2.k  (437> 


Equation  (437)  must  be  satisfied  for  any  values  of  r  ,  ?p  ,  and  s  . 
Functions  separated  In  this  manner  on  the  two  sides  are  to  be  separately 
equal.  Equating  parts  In  a  and  cp  , 


Znk(,)  ’  Z2nk  *  ,lnk*  +  *onk) 
*nk(<p)  *  »2nk  '  ^K^onk) 


(438) 


(439) 
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giving 


M 

z_,(«> 


2nk<*> 


2 


and 


0($) 


“  (6k) 


Than  equating  the  parts  in  r  on  both  of  tha  aids*  of  aquation  (430), 


^  Enk 


+  *"  +j*.isiu!r 

nk  r  2  *nk 


R  U 

-at 

2 

r 


-  R 


2nk 


(440) 


Subatituting  for  froai  equation  (430),  aquation  (440)  becomes 


•  h  • 

,  R  .  I  L 

‘ 

Li  M  r 
^n  n 


♦  R  li  +  MTt  i 

nk  r  2  2 

^  r 


*  -^[w6kKr)+»6„kHtt(‘v)]  «*» 


Let 


C**r.) 


Then  aquation  (441)  takas  tha  following  form. 


Uiat  +  1  Hit  +  (i  .  I  , 

dr2  r  dr  r,  I  nk 
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(443) 


*  ^  [*7  ♦  *,  Cw] 


where 


12  k  B5„v 

B«  •  —  ■  -  I?  iu  ■ 

7  i  12  u  B5nfc 


(444) 


and 


12  k  B.  . 

b8  -  — r *  - 12  tk  »*„„ 


(445) 


Let 


Rnk  “  R(Tl)  “  U(Tl)  J  / - W  *  V(t9  H(l> _ (tt)  (446) 

V  36k2  +  1  V 36k2  ♦  1 


be  the  particular  Integral  of  equation  (443). 
Then 


R'OT)  -  U '  (r))  J  _ (n)  ♦  UfrO  J1 

V  36k2  +  1  V36k2  ♦  1 


V'(tt)  H(l2 


(1>\ 


Y - 0l>  ♦  V(tt)  (n>  (447) 

▼  36k2  ♦  1 


Take 


U'(ri)  J  ^ - +  V'  (r|)  H - (n)  -  0  (448) 

V 36k2  +  1  V 36k2  ♦  1 
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Then  differentiating  equation  (447), 


I"(tt)  -  U’(ri)  J’  _ (rO  ♦  U(n)  r  - (tt)  + 

V36k2  ♦  1  V 36k2  ♦  1 


V'(ri)  H(l*! _ (rtf  +  V(rO  H(l^] _ _ (rj)  (449) 


V36k2  M 


Vttk2  ♦  1 


Using  equations  (446)  through  (449)  in  equation  (443)  and  rearranging 

gives 


«w  |j". - w  v. - 

L  Vj*kJ  +  i  Vatk2  +  i 


Or)  ♦ 


1  «<!)' 


U'Ol)  J' 


1  J  r- 

-  (n)  ♦  1 

W  H'*' 

1  »3*kz  +  1  J 

L  yyt* 

(rj>  ♦ 

L .  at2  .t  il 

[  ^  J 

|h<1> _ 

'  V36k2  ♦ 

_ (ti) 

♦  V'  (n)  H(ll' 

(n) 

♦  1 

Vj6k2  ♦  l 

-  A  Fb,  , 

.(ri)  ♦ 


W* + ' 


» <?«>] 


(450) 
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•r 


Noting  that  J  _ _ and  (r)  art  solutions  of  ths 

V 36k2  ♦  1 
Bassal  aquation, 

y”(t i)  +  I  y  .  o 

rssults  in 


rri—™  +  n  yn— w  + 

V36k2  ♦  1  Y36k2  ♦  1 


1 


36k2  ♦  1 
2 


r\ 


J  , - (r)  • 

Y36k2  ♦  1 


and 


— (*n)  ♦ 
♦  1 


1 

r\ 


« <i>_ 
^Mc2 


— Or)  ♦ 

*  l 


i 


Using  aquations  (452)  and  (453)  in  aquation  (454), 

»' (n>  r  - tn>  +v’(r1)  h(1>! _ <r,) 

V36k2  +  1  V36k2  +  j 

‘  ■jf#7J6kW  +  BeH»i>W 
115 


following 


(451) 


0  (452) 


0  (453) 


(454) 


Solving  for  U'  (tj)  and  V’ (r)  from  equations  (448)  and  (454), 


bOCAUM 

Vn)  *pUw  •  yi>  n‘l)'<ti)  •  ^  (*”> 


Therefore,  aquations  (455)  and  (458)  bacoaa,  raspactlvaly, 


U'(ri) 


¥  I- 


♦  » 


(458) 


v’(n> 


■  *  ^{'v£?71w[‘’  J"w  *  ■■  c'i 


(459) 


Integrating  the  last  two  expression*  with  respect  to  r)  , 


U  *  T  |*7  *1  +  *8  x2) 


and 


V  *  -  T  1*7  *3  4  »8  h) 


where  I, 


dri 


dr| 


(460) 


(461) 


(462) 


(463) 


(464) 


(465) 


Equations  (462)  through  (465)  aay  be  represented  In  general  by  the 
following: 


I 


(466) 
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vhoro 


Z  and  1  stand  for  or  J 


Y  stand  for  V36k2  +1 
0  stand  for  6k 


Z^  stand  for  Z^(r)> 


Zp  stand  for  Z^Ol) 


> 


(467) 


Mow 


fa  .  Vi  *  Vi 

n  2  l 


(468) 


Zhan  aquation  (466)  bacoass 


Tho  following  intagratlon  formula  Is  to  bo  notad: 


/i 


.  z_  .  7  -  z  T  .  z  T 

Ail  dt,  -  „ .  _a_a 

*1  P  4  »2-«2  P  ♦  4 


p  -  q 


(469) 


(470) 


Than 


J$\V  i4”  ’ 


z  .  T.  .  -  z  Y.  ,  z  F.  . 

y-l  8-1 - v  8-2  .  V  8-1 

^  2  .  v2  y  ♦ 

Y  -  (|  -  l)  T  w 


(471) 


118 


1 


md 


!< 


i  *  i  t...  -  z  r.  z  r 

i  *v  Vl  4t>  •  Y:2  - V1  * 

r)  v  W  Y*  -  (p  +  l)2 


11  y  fr**  1  "  Y2.(p+l)2  Y+P+l 

Using  aquations  (471)  And  (472)  In  AquAdon  (469)  glvas 


2  ,  2.  ,  *  Z  Z.  „ 

v-1  6-1  v  6-2 

\-lVl  *  2y7.\ 

M 

m 

5  i  “  ■  *• 

i  Y  •  •  1) 

Y2  -  (»  +  l)2  1 

IzL 


Z_. 

♦  — 


Y  vn-l  Y  +  IM 


(473) 


Froa  Aquation  (467),  substituting  for  y  6  * 


Y2  -  (6  -  l)2  -  (^6k2  +  l)2  -  (6k  -  l)2  •  12  k 
Y2  -  (0  ♦  l)2  -  (V^k2  ♦  ij  -  (6k  -  l)2  -  -  12 


Than  tha  foliowing  sra  aaslly  obtAlnad: 


:  -i  -»-i  •  zv  z.-2 1  vi  m  ~  s ;« 

Y  *  tt  *  »2  Y2  *  (»  +  l)2 


-  z.  z . 


.  Vi  ( Vi  '  Vil  *  *y(  *»  •  ri-a1 


12  k 


(474) 
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Using  aquation  (474)  In  aquation  (473),  than  substituting  for  y  and 
0  tram  aquation  (447)  raaulta  in  tha  following  aquation: 


Equation  (473)  is  aquivalant  to  aquation  (466)  with  its  integration  accoa- 
pllshad.  By  analogy,  similar  axprssslons  aquivalant  to  aquation  (462) 
through  (463)  ara  as  follows  (in  sequence). 
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Tha  total  solution  of  la 

*i  ■  *1  <»•*•»>  *  ?i  <*»°) 

Equation  (480)  contains  g^(r,9,a)  ,  tha  solution  of  tha  homogeneous  part 
of  aquation  (408)  given  In  aquation  (427).  Using  aquation  (438),  tha 
following  la  restated: 

*.kw  *■*<»>  «*•»> 

Znk(a)  and  #nk(f)  *r*  given  In  equations  (438)  and  (439),  and  restatad 
herein: 

*»k(,)  *  ***  ■  •uk,  +  ,od 
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«nk(<P)  ■  «2»k  *  H4k*+%J 

The  most  difficult  value  to  establish  is  R(r)  ,  ss  shown  in  the 
analysis  contained  in  equations  (440)  through  (479).  The  subsequent 
steps  nust  be  followed  to  obtain  R(r)  :  (1)  calculate  first  the 
Bessel  expression  in  equations  (476)  through  (479);  (2)  Introduce  the 
obtained  to  Into  equations  (460)  and  (461)  having,  then, 

U  and  F  .  Entering  this  expression  Into  equation  (446),  gives  the 
values  for  Rftk(r)  : 

a  -  T(n)  -  u(n)  j  __(n)  +  vot)  h  — . - jr) 

11  V36k2  +  1  Y36k2  *►  1 

Having  the  values  of  distortion  ^  »  then,  the  stress  distribution 

Is  deteralned  by  equations  (1)  through  (5)  and  equation  (7)  for  a  . 

Thus,  the  aatheaatlcal  -  physical  fundaaantals  for  the  Internal  aachanlcs 
of  parallel  fibers  are  established. 

Review  of  equations  (436)  and  (439)  clearly  indicates  that  the  ('•  are 
of  an  undulatory  nature.  But  since  the  solutions  are  suns  over  several 
eigenvalues  characterised  by  paraaaters  n  ,  they  could  be  regarded  as 
a  Fourier  series.  The  coefficients  of  the  series  obviously  show  a  certain 
characteristic  sbout  a  value  nQ  .  The  assus^tlon  can  therefore  be  that 
the  functions  g  consist  prlasrlly  of  those  cosine  functions  whose  order 
la  grouped  within  a  narrow  region  around  this  value  nQ  .  The  conclusion, 
then,  can  be  that  the  functions  {  will  consist,  for  the  aost  part,  of 
the  functions 


a  ns 
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A  b 

But  since  n  Is  s  function  of  r*  end  ■£>  end  the  aster  lei  chare c- 
o  L  b 

o 

terlstlcs,  the  wavs  length  will  elso  depend  on  those  vs lues,  besides 
showing  special  dependence  on  the  fiber  disaster. 

A  photoalcrogreph  of  residual  stress  In  Epon  820  resin  surrounding  a 
saaller  disaster  E-glass  single  fiber  Is  shown  In  Figure  4.  Megnlflcatlon 
Is  100X.  Fiber  disaster  was  0.0004  In.  This  free  fiber  was  pieced  In  e 
drop  of  resin  end  first  cured  et  250*F  for  30  sdnutes.  Examination  et 
75*F  In  polarised  light  following  this  cure  showed  no  residual  stress, 
end  the  fiber  was  straight.  However,  after  an  additional  cure  of 
30  alnutee  et  350*F,  the  seaple  appeared  as  shown. 


Figure  4.  Photoalcrogreph  of  a  Fine  Fiber  In  Resin 
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APPENDIX 


DERIVATION  OF  FUNDAMENTAL  EQUATIONS 

The  following  la  a  brief  discussion  of  the  basic  physical  relations 
which  will  be  required  in  later  applications  to  a  fiber  reinforced  matrix 
under  various  loads.  In  view  of  confusion  in  the  literature  about  whet 
is  and  what  is  not  a  tensor  or  tensor  component,  this  introduction  is 
considered  necessary  to  define  terminology. 

First  Order  Theory  Introduction 

Consider  the  transformation  from  an  orthogonal  Cartesian  coordinate 
system  to  a  mere  general  curvilinear  system  of  coordinates.  The  trans¬ 
formation  is  a  one-to-one  transformation  all  of  whose  derivatives  exist: 


^  =  (« 

Let  be  a  set  of  orthogonal  unit  vectors  so  that  the  scalar  product 

gives 


X  •  T; 


^ J  9 

J  « 1  ,  -  J 


(2) 


6.  . is  defined  by  the  scalar  triple  product  (Reference  1,  pages  16 

ani  20) 


T.  x  I;  -  lie 

x-  J 


(3) 


so  that  the  vector  product  is  given  by 


<*) 


^/Lj  n  ■  0,  for  any  equal  subscripts 

■  1,  for  cyclic  interchange  of  subscripts  1,  2,  3 

■  -1,  for  non-cycllc  Interchange  of  subscripts 


The  upper  index  jX  used  to  conform  to  subsequent  tensor  notation. 
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Consider  Che  vector 


Cartesian  coordinates 


whose  differential  displacement  is 


^5  -  AT  r  I, 

> 

repeated  indices  sunmed  which  becomes  under  the  transformation  (1) 


d 


X 


> 


(5) 


(6) 


where 


> 


(7) 


Equations 

formations 

tangent  to 


(S)  and  (7)  represent  the  contravariant  and  covariant  trrns- 
respectively.  The  vectors  are  the  metric  vectors  tangent 

dJ. 


the  curvilinear  coordinate  elements 


The  differential  arc  length  da  la  given  by  the  scalar  product 


ds*  *  dR  -  df  *  (/*-*•  J»j  T  ■  J. 

-  dx*  dxj  a  •  . 

-  Ju*  duj 

=  Ju^du^  *  3j  ■*  du"  du-  • 

^  J 

where  g^  ,  the  metric  tensor  la  defined  by 


(8) 


3 


=  r  StL  ^4 


(9) 


Although  there  le  no  numerical  distinction  between  coverlent  tensor 
end  contreverlent  tenaor  quantities  In  Cartesian  coordinates,  the  mainte¬ 
nance  of  the  upper  end  lower  Indices  Is  a  convenient  accounting  arrange¬ 
ment.  For  example,  If  the  Indices  In  equation  (7)  ere  raised,  the 
contreverlent  metric  vector  system  is 


I 


From  these,  one  may  define 


(7a) 


3 


(10) 
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and  show  by  aquations  (7),  (7a),  and  (10)  that 

jk  '  3j  »  £j 

3ki3:i ;  *  6j 


(10a) 


By  vlrtua  of  thasa  aquations  and  othar  propartlas  of  tha  natrlc 

•  i 

tansor  componants  •  thasa  may  ba  utilisad  for  raising  and 

lowering  indices.  ^ 


A'.S^Aj  ;  Aj-Su* 


Note  from  aquation  (9)  that  tha  absolute  value  of 


and  hanca 


3el 


Si 


are  unit  vectors. 


Dlffarantlstlng  aquation  (9), 


J4-  . 


(>£,  ^  j 

_  t  '  3k  +  3j  *  „  l 


du  ou 


»  & 


ja!al 


-  +  G  ^ 

•  *  r  S 


Jtm  _  U. 


du  du 


Interchanging  i  ,  j  ,  /c  ,  adding  and  subtracting,  yields 


(14) 


J 


which  It  the  Crlstoffel  symbol  of  the  first  kind.  The  Crlstoffsl  symbol 
of  ths  second  kind  Is  doflnsd  by 


{d}  -  / 


Also, 

ar 


du 


■*-*-  r 


dukdu^  *  da*  Da 


2^4  T<"  £ 


(15s) 


£ 

du^du-* 

*  •  9 

aaj 

*  Ck  \  j]  : 

-  £*}  Sj 

The  left-hand  side  of  this  equation  Is 
covarlant  derivative: 

*.•5-1 

defined,  i 

(15b) 


(16) 


In  similar  manner,  the  following  can  be  shown: 


3^  - 


(16s) 


end  from  equation  (14), 
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*5 

> k  du 


?  '  ft)  3<j '  ft} Sil  =  0 


(17) 


in  similar  manner, 


3kJ 


(17a) 


ao  chat  cha  covariant  derivations  of  the  metric  vectors  and  the  metric 
tensor  are  aero. 


The  determinant  of  from  equation  (9),  becomes 


3u~  da 


4 

du* 


(18) 


X 


(19) 


3  *  l|3ikll  * 

and  by  equation  (7) 

3  *  Ofr  *  £  *  5jJ 

Hence,  by  equation  (3), 

5r  •  3 ’  3*.  =  Jr  *  3s  '  3*  -  erJ>i  Vg"  (20) 

1  Si  ^ 

Fran  these  and  the  fact  that  ^  #  <Jj  *  u  J  > 

3r  ' 


3* 


(21) 


and  likewise 


where 


rst 

F 

>~wt 


9-t 

and  fc***  -  e 


(22) 


rvt  ^ 


are  tensors. 

The  area  of  the  triangle  determined  by  two  of  the  three  vectors 


"  5t  <L  d  (^J)  not  summed 


(23) 
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if 


dsi 


u  * 
J 


*  i  fe  *  *j]  ^  A 

=  i  [?k]^r 

s  X  ^  dlu.* 


(24) 


By  equation  (21), 

?i  '£& 

la  a  unit  vac  tor  normal  to  tha  araal  clamant  fll$^  ;  hanca, 

*dS  •  f  ife  *Si 


(25) 


where  tha  aum  of  tha  vector  araaa  of  tha  facea  of  tha  tetrahedron,  formed 
by  tha  metric  vectora  d 5^  ia  equal  to  tha  vector  area  of  tha  fourth 

face "nd  5  ,  where  n  la  tha  unit  normal  to  tha  fouxth  face  d S  . 

Batting 


n  -n,  3 


Ut) 


one  obtalna 


(26a) 
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Fir  it  OrdT  Strain 


Consider  a  displacement  vector 


v  •  v‘£  =  v*  fsZZ  •  ~7  -  r&x 

3^ 


(27) 


where  the  ere  the  physical  components  of  the  displacement  vector 


^  ,  since  ere  now  unit  vectors.  are  not  components  of  a 

tensor  and  are  thus  confined  to  a  particular  coordinate  system.  Taking 
the  differential  of  equation  (27) 


i~  ■  1*4  I,  ♦  a 


\Da 


duJ 


i  Ari 


(28) 


which  by  virtue  of  equation  (ISb)  becomes 


■( 


du  J 


and  by  equation  (16) 


*  vj 

*  V-J  IA 


Because  V  •  • 

J 


W*>J 

is  not  syometrlc,  the  following  may  be  written: 


X>J 


iul  +  z*L± 

z  z 


(29) 


z  D:  ;  +  <^± 


so  that  the  total  displacement  becomes 


"5  *  V  +  J*  du' 

J  J 


(30) 
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Here ,  the  first  and  third  terms  represent  s  pure  displacement  and  a  pure 

rotation  respectively  and  the  second  term  IX »  ^  du **  represents  the 

strain  dlsplacesMnt  where  P-  '  is  the  strain  tensor.  The  rotetlon 

J 

tensor  Q:  ;  is  not  considered  here,  but  may  appear  later  where  local 
J 

couples  give  rise  to  nonsero  values. 


The  Stress  Tensor 

Consider  an  elemental  tetrahedron  enclosed  by  the  coordinate 
surfaces  dS  :  snd  an  arsa  dS  *  which  is  in  equilibrium  under 

forces  acting  across  these  faces;  i.e., 

JdS  *  2  dS-  (31) 

w  *  A* 

4* 

where  ^  is  the  stress  vector.  ^  is  an  invariant,  but  “  are 
not,  in  the  sense  that  they  act  across  a  particular  set  of  coordinate 
face  elementa.  Prom  equations  (26a)  and  (31), 


n.  O',  yq  .7  (32) 

Ji  x. 

By  the  definition  of  * .  equation  (25),  these  are  covariant  components 

and  hence,  because  of  the  lnvsrlance  of  O'  .  the  following  may  be 
written: 


CT  { 3  «  ( T  (33) 

•  t 

which  are  now  contravar lent  components  of  a  tensor,  and  O'  J  L»  defined 
as  the  contravarlant  of  the  stress  tensor. 

Prom  equations  (32)  and  (33) 

?  -  dr* J  pu  (34) 

4lso, 

*5  *  <rJgj  os) 
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so  that 


(7  *  <T4'Jn^=3J*cr  06) 

where  CT  ^  arc  the  contravariant  vector  compcnenta  of  the  stren  normal 
to  d  Sj  . 

The  three  atreia  vectora  may  be  written  (Reference  1,  p.  78) 


(37) 


To j  <»»> 

*  I 

where  define  the  physical  components  of  the  stress  tensor „  but 

refer  to  a  specific  coordinate  system  and  are  not  components  of  a  tensor. 


Equilibrium  Relations 


In  Figure  3,  the  stress  vectors  acting  across  the  face  of  the 
elcmenta  of  the  parallelepiped  formed  by  the  three  vectors  q  .  du 
are,  referring  to  Figure  3, 


<±) 


-or.  dS  ; 

-o  (-»■)  I 


or 


& 

-c  (i) 


A» 
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H 


If  no  couples  set  on  this  element  of  volume,  then  from  equation  (39), 
keeping  only  first  order  terms, 


(40) 


or  from  equations  (37)  end  (22), 


e. 


k 


<r  x  j 


0 


from  which  it  follows  thst 


•  •  '  ; 

*  crJ 


(40a) 


(41) 


and  the  stress  tensor  is  symmetric. 


The  remaining  equilibrium  equation  states  that  the  sum  of  the  forces 
acting  on  the  element  of  volume  are  aero.  Thus,  if  the  mass  accelera¬ 
tion  of  the  volusir  is  -  ?  f  d  V  «nd  the  body  forces  are  PrdV  , 

where  9 


*  [91  *  32  9j]  d*1  Ju 1  Ju 


z  if  Ju1  da*  clu3 


(42) 


then  the  equilibrium  of  forces  is  expressed  by 


0 


or  by  equation  (37) 

* 


33 U 


C>0L 


+  PFJ*T  3;  *  pf 


which  can  be  written 


(43) 


(44) 
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but 


C )U 


♦if? 


(44a) 


1  3/j~  f  r  1 

*3  0^  =  M 


and  from  equation  (15b), 


£♦£} 


do*^ 
a 


a.r 


fc) 


Lj  f.  ppJ  *  pf J 


or 


o;N  tpF‘» 

J 

•  I 

Since  Q'^J  la  a  symmetric  tensor,  the  contracted  covariant 
derivative  may  be  written  in  the  form 


<r  • 

j 


At* 


(45) 


(46) 


(47) 


which  is  much  simpler  to  evaluate  for  a  particular  coordinate  system. 


Stress-Strain  Relations 

The  conditions  of  equilibrium  (41)  and  (45)  furnish  six  relations 
for  tl '  determination  of  the  nine  stress  tensor  components.  The  three 
remaining  are  furnished  through  the  physical  equation  functionally 
relating  the  stress  and  strain  tensors.  This  functional  relationship 
may  be  expressed 

art)  M.  <*»> 

where  the  functions  are  expendable  in  a  Taylor's  series.  If  this  is 
expanded,  retaining  only  first  order  term*,  the  result  is  a  statement 
of  Hooke's  law. 
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(49) 


<rxJ  -  +  Pi  Dri 

or  Inversely 


*  cn  =  cn^j 


>  (50) 


where  the  constant  tarma  Cpe^  >  A  may  rapraaant  Initial  stress- 

strain  conditions  due  to  polymeric  shrinkage  and  differential  thermal 
expansion. 


The  sysseetry  properties  of  a  material  are  independent  of  coordinate 
choice,  so  that  from  the  evaluation  of  the  Hooke's  lev  elastic  tensor  a 
Cartesian  freras  of  reference  may  be  chosen. 


In  isotropic  media,  such  as  the  reinforcement  or  the  matrix,  the 
elastic  tensor  P%  is  Invariant  under  rotation. 

A^n  „  kr  A;  ^  A^kl 


where  A  \  la  a  transformation  of  pure  rotation  and  hence 
J 


a!  A} 

~  J 


G 


Cr  •  • 


It  follows  that 

•  O 

+  rj&*x 


(51) 


(52) 


(53) 
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where  X 
of  &  -*j 


fA  j  ?J  ere  eceler  multipliers, 
end  P  . 


In  view  of  the  symmetry 


A 


(54) 


which  leeds  to 


M  M 

end  hence 

/\^kl  =  \G^Gkl  +  M 


(55) 


so  thet  equetlon  (49)  becomes,  piecing  =  0  where  X  f  fiA.  are 

Leme's  constents 


•  f 

cr 


&3*J3ki  +  *  +6*SGjk)]pk£  (56> 


In  generelised  coordlnetes, 

G"*j  =  [xgV*  +  ^  (a^3jl  +5X<3Jk)j  D*i 


(57) 


Plec lng 


E  v1 


) 


A< 


where  E  is  Young's  modulus  of  elestlclty  end  V  1*  Poisson's  retlo. 


The  following  is  obtained  (Reference  1,  p.  162): 


1 

tr‘J*  f  3'V"  » 

<r*J*-rf5-  CDiJ  3*J'pi) 

°r  *\ 

*  lf7  (D-j  +  WV  9ij  pi) 

0-4  s  -1-  (pi  +  — i!~  <$-f  d}) 

°j  i*v  \  J  i-2>>  J  ll  J 


Inverting, 


D^j  *  i±Jd  O'^j  -  ^  <W 


-  ii'Dimiw  rrm  EMMZimn  ttti 


Tha  conatant  tarma  In  aquation  (49a)  can  ba  utilised  to  rapraaant 
the  Internal  affects  of  polymerisation  shrinkage  and  thermal  expansion. 

If  0 C  repreaents  the  linear  coefficient  of  thermal  expansion  and 
the  polymeric  shrinkage,  then 


cj  =  3j  («T  *A) 
*  <5j  (*r  +>3) 


I 


and 


D?  -  h\  (ofT  +  >3)  * 

J  J  fc 


or  Inverting 

0 4  =  JL1 

pf  +  df  | 

J  J 

to 

f  V 

iwv 

o-Aj  *  JL . 

f* 

*  •  *  • 

Lp4j<-3*jl 

[hv 

Pf  -  I® 

*  1 


( 


Differential  Equations  for  the  Displacement  Vector 

Equation  (46),  for  static  applied  forces,  become* 


(T^J  4-  fF  *0  (64) 

0 

Prom  equation  (63), 


Contracting  O'4'-*  to  (y 

e— 
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e-ear 


i 


(66) 


Equation  (64)  now  becomes 


ypj  I!*  .  0  (67) 


or  lowering  indices 


(68) 


where  F^  =■  F^  can  *P  +  Fz  vo 


F}  f*  S»n  ‘f  f 


fj 

Fx 


are  body  forces 


r  cos  <p 


in  Cartesian  coordinates 


whereas  F. 

A* 


are  in  cylindrical  coordinates. 


146 


u*  cos  ux 


u1  5<A  u*  }  UX  • 


Than 


3-j  *  0  , 


3n  * 

1 

•  9*z  *  («*)  ,  3SS  *  1 

3*  * 

and 

1 

'  3ai  *  r*  >  JiJ  *  ^ 

#  t 

3*J  3 

0  , 

^  ^  j 

3*  3 

1  , 

„  _  1  Ji _  , 

9  pp  >  3  -1 

3 11  m 

i  , 

32i  *  ^  '  9***1 

Tha  only  nonaaro  Criatoffal  synbola  of  tha  aacond  kind 

{A}  *  - ' 

{*»}  3  {«}  3  r 

(70a) 


I 


Further,  by  equation  (49) 


where  by  equation  (27) 


Hence: 


D 


D 


DIS  - 
b«  = 

Di5  = 
D,J  = 


Jdl 

3* 

2  (_r  d<p  ar  r  ^ij 

X  fiSf  +  1  iSfl 

*  1.3*  ^  39  J 

x  [a-  +  -as.1] 

a  L»r  3t J 
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Utilising  equations  (27)  and  (47)  in  equation  (68) 


-i-K  +.  _i_  +  -1- i  *!*t-  + 

dr1  <Jr*  8?*  i  dl7  T  2(l-2»)  r  drdt  r 

j  _  j»f »  1  . 

vd*  *(*-*>0  7*  d<f 


E 


C*S<P  - 


Sjl  *'<*■<()  , 


# 


* 


x  ifk 

&r*  ^  1-iV  r*  8<f* 


i  ,  i.  i!$i.  +  _J — .  J. 
T7F2Z)  r  Drtl  ^  2( i-M  r 


d*tdr 


( 


*r*  i  "j  v 


(Si  (ini’  «■  3*  c«4  «  0 


0  (73) 


(74) 
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1  ±!ia.  +  JL  +  JLz )L  4  i 

2  2c*  2c»  Vf*  T  Wv>  >**  +  2(l~2v)  "a^I 


_  ,  -  j  -  l*L 

l(i-2v)  ZVdz  2(1  -*w  r  a* 


+  -L 

2r 


j*a 

dr 


JUL 

1-2W 


('S‘gN 


jLn> 

c 


(75) 
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